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Abstract 


The  ship  motion  computer  program  SHIPMO  has  been  revised  to  increase  its  capabilities  in 
several  different  areas.  The  new  version  of  the  program  can  compute  forces  relative  to  local  shij) 
axes  and  the  resulting  incidence  of  motion-induced  interruptions  or  sliding  events.  Seakceping 
positions  may  now  be  offset  from  the  ship  centreline.  The  program  now  includes  expanded 
wave  spectra  capabilities,  including  multi-directional  spectra.  Motion  computations  may  include 
hydrodynamic  end-effect  terms,  which  can  be  important  for  transom  stern  ships  at  forward 
speeds.  SHIPM05  also  has  an  option  for  predicting  propeller  emergence,  which  can  influence 
added  power  requirements. 


Resume 


SHIPMO,  le  logiciel  de  simulation  des  mouvements  d’un  navire,  a  etc  revise  de  maniere 
a  accroitre  ses  possibilites  sur  divers  plans.  La  nouvelle  version  du  logiciel  peut  calculcr  los 
forces  par  rapport  aux  axes  locaux  du  navire,  ainsi  que  la  frequence  resultante  des  evenements 
de  glissement  ou  des  interruptions  dues  au  mouvement.  Les  positions  de  tenue  en  mer  jreuvent 
maintenant  etre  decalees  par  rapport  a  Taxe  longitudinal  du  navire.  Le  programme  comporte 
maintenant  des  fonctions  et endues  de  calcul  des  spectres  d’ondes,  y  compris  les  spectres  muhi- 
directionnels.  Les  calculs  du  mouvement  peuvent  tenir  compte  des  effets  hydrodynamiquc^ 
finals,  qui  sont  parfois  importants  dans  le  cas  des  navires  a  tableau  arriere  en  marcho  avant. 
SHIPM05  comporte  aussi  une  option  de  prediction  de  I’emergence  de  I’helice.  un  evencinent  qui 
peut  accroitre  la  consommation  d’energie  du  navire. 
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Notation 


A{P,) 

A, 

A44 

a 

ai»02»a3 

B 

Cb 

CG 

F, 

Flat 

F long 

F, 

fr 

GM 

9 

H, 

h 

^U' 

fCxi 

L 

LCG 

1{F 

M.M 

Mp 

Mil 


m„ 

n, 

0^ 

F 

^port 

^stbd 

5(T).5(a.v.) 

S{u;^,,u) 

Ssiu^) 


T 

Taft 

Tfort 


Tport 

Tt 

Tr 

F 


spectral  normalization  factor 
sectional  area 
added  moment  of  inertia 
wave  amplitude 

accelerations  in  global  axes  at  position  on  ship 
beam 

block  coefficient 

centre  of  gravity 

excitation  force  for  mode  i 

lateral  force  estimator 

longitudinal  force  estimator 

vertical  force  estimator 

sectional  excitation  force  for  mode  i 

metacentric  height 

gravitational  acceleration 

significant  wave  height 

height  above  deck  of  CG  of  person 

wave  number 

roll  radius  of  gyration 

ship  length  between  perpendiculars 

longitudinal  centre  of  gravity 

horizontal  distance  from  CG  of  person  to  tipping  point 
spectral  directional  spreading  function 
number  of  slides  per  second  in  port  direction 
number  of  motion-induced  interruptions  during  task 
spectral  moment 

component  of  normal  pointing  outward  from  hull 
tipping  point 

spectral  directional  spreading  parameter 

port  sliding  force  estimator 

starboard  sliding  force  estimator 

unidirectional  wave  spectral  energy  density 

directional  wave  spectral  energy  density 

spectral  density  of  sliding  estimator 

half  stance  width  of  person  on  deck 

draft  or  wave  period 

aft  tipping  estimator 

forward  tipping  estimator 

peak  wave  period 

port  tipping  estimator 

task  duration 

sectional  draft 

period  corresponding  to  energ,v  averaged  wave  frequency 


T-\  energy  averaged  wave  period 

Vk  ship  speed  in  knots 

V  effective  flow  velocity 

I,  y,  z  coordinate  system 

a  exponent  for  JONSVVAP  spectrum 

0a  direction  of  wave  propagation  relative  to  ship  velocity;  0^  =  180°  for  head  seas 

X  compass  bearing  of  ship  forward  velocity 

6  azimuth  angle  of  tipping  point  0$ 

ship  angular  displacements 
r()  gamma  function 

A;  spectral  shape  parameter 

A,  tank  valve  resistance  coefficient 

fis  static  coeflficient  of  friction 

v  compass  bearing  from  which  wave  component  is  approaching 

I'mi  mean  compass  direction  (from)  of  spectral  energy  for  component  i 

p  water  density 

<7  JONSVVAP  spectrum  coefficient 

w  wave  encounter  frequency 

Wp  peak  wave  frequency 

Wu,  wave  frequency 

w_i  frequency  corresponding  to  energy  averaged  wave  period 

Q  spectral  amplitude  parameter 

A  ship  mass 


VI 


1  Introduction 


The  ship  motion  computer  program  SHIPMO  (1,  2,  3,  4]  has  been  revised  to  increase  its 
capabilities  in  ship  seakeeping  analysis.  The  changes  have  been  implemented  in  a  number  of 
different  aspects  of  the  program. 

To  provide  improved  estimates  of  ship  operability,  SHIPM05  can  predict  vertical,  lateral, 
and  longitudinal  motions  at  off-centre  seakeeping  positions.  These  motions  can  subsequently 
be  used  to  compute  lateral  and  longitudinal  force  estimators  and  the  resulting  incidence  of 
motion-induced  interruptions  (Mil’s)  or  sliding  events  at  seakeeping  positions  using  the  theory 
described  by  Graham  [5].  Computed  lateral  and  longitudinal  force  estimators  which  result  from 
accelerations  relative  to  the  local  ship  axes  can  also  be  used  in  ship  systems  design  as  described 
by  Beck  et  al.  [6]. 

Several  new  wave  spectra  have  been  added  to  SHIPMO.  The  JONSWAP  spectrum  [7]  is 
suitable  for  fetch-limited  conditions.  The  Ochi  and  Hubble  6  parameter  spectrum  [8]  is  a  uni¬ 
directional  spectrum  having  separate  swell  and  sea  components.  The  10  parameter  spectrum 
[9,  10,  11],  which  was  originally  proposed  by  Hogbcn  and  Cobb  [12]  as  a  multi-directional 
extension  of  the  Ochi  and  Hubble  6  parameter  spectrum,  is  also  available.  The  new  program 
can  also  use  multi-directional  hindcast  spectra  given  by  the  user  or  from  the  Offshore  Data 
Gathering  Program  (ODGP)  hindcast  model  [9,  10].  Hydrodynamic  end-effect  terms  [13]  on 
transom  sterns  are  available  as  optional  computations.  Incorporating  end-effect  terms  gives 
improved  motion  predictions  in  many  cases  for  transom  stern  ships  at  high  forward  speed  [14]. 
A  revised  method  for  computing  surge  excitation  should  give  improved  force  predictioits  along 
the  hull  and  includes  forces  on  the  transom.  S1IIPM05  can  also  predict  propeller  emergence 
frequency,  which  can  influence  added  power  requirements. 

To  simplify  program  input  for  the  new  version  of  SHIPMO,  vertical  coordinates  of  the 
centre  of  gravity,  deck  edge,  anti-rolling  tank,  and  scakeeping  positions  are  input  relative  to  the 
baseline;  however,  input  files  for  older  versions  do  not  need  to  be  modified  to  run  with  version 
5  because  the  new  version  is  capable  of  recognizing  earlier  version  files  and  making  the  required 
adjustments  to  the  input  data. 

SHIPM05,  like  SH1PM04  [4],  does  not  have  a  wave  load  capability;  thus.  SH1PM02  [2] 
must  be  used  if  wave  loads  are  required. 


2  General  Discussion 

This  section  discusses  several  important  areas  to  note  before  running  SHIPMO.  The  user 
may  also  wish  to  consult  the  theoretical  basis  of  the  program,  given  in  Refeiences  1.  3.  5,  and 
13. 

2.1  Coordinate  Systems  for  Ship  and  Waves 


The  coordinate  system  used  for  the  motion  calculations  is  illustrated  in  Figure  1.  Its  origin 
is  located  at  the  equilibrium  position  of  the  ship  centre  of  gravity  (CG).  The  x  axis  points  to 


the  bow  and  y  axis  extends  to  port,  and  both  axes  are  parallel  to  the  undisturbed  free  surface. 
The  z  axis  points  vertically  upward.  The  coordinate  system  translates  with  the  mean  velocity 
of  the  ship,  maintaining  a  fixed  orientation  with  respect  to  the  free  surface. 

The  sea  direction  /?,  is  defined  in  Figure  2  as  the  angle  between  the  mean  velocity  vector 
of  the  ship  and  the  direction  of  propagation  of  the  wave  train,  which  is  normal  to  the  crests  of 
the  waves.  Thus,  /J,  =  0°  for  following  seas,  fia  =  90°  for  beam  seas  from  port,  and  [is  =  180° 
for  head  seas. 

For  the  10  parameter  and  hindcast  directional  spectra,  an  energy  component  is  dependent 
on  the  compass  bearing  v  from  which  the  component  is  approaching  as  defined  in  Figure  3.  If 
the  ship  forw'ard  velocity  has  a  compass  bearing  of  Xi  then  the  relative  sea  direction  Z?*  of  the 
energy  component  can  be  determined  as  follows; 

Us  =  I/  +  180°  -  X  (2.1) 

For  the  input  hindcast  spectrum,  the  energy  in  each  directional  bin  is  the  wave  energy 
going  to  that  direction.  This  convention  is  changed  within  SHIPMO  such  that  the  energy  in 
each  directional  bin  is  the  wave  energy  coming  from  that  direction. 

2.2  Accuracy  of  Hull  Form  Definition 

Hull  form  input  for  SHIPMO  can  be  in  the  form  of  sectional  offsets,  or  sectional  beam, 
draft  and  area  coefficient.  Good  roll  predictions  require  accurate  offset  data,  since  roll  damping 
is  very  sensitive  to  section  geometry.  Therefore,  it  is  recommended  that  the  user  use  offsets  to 
describe  the  section  geometry  whenever  the  lateral  plane  responses  are  required. 

Sectional  hydrodynamic  properties  can  be  computed  using  either  the  conformal  mapping 
method  or  the  close-fit  method.  It  is  possible  to  input  offsets  and  still  use  the  conformal  mapping 
method  to  take  advantage  of  faster  computation  times.  When  the  conformal  method  is  used 
outside  the  range  indicated  in  Figure  4,  which  has  been  adapted  from  Reference  1-5.  an  error 
message  will  indicate  which  stations  are  outside  the  appropriate  range.  In  general,  the  close-fit 
method  is  recommended  whenever  this  is  the  case;  however,  if  only  bow  stations  having  very 
small  area  coefficients  (sucli  as  those  occurring  on  some  warships)  are  involved,  the  added  mass 
of  a  station  will  be  a  small  fraction  of  the  total  added  mass  and  the  error  introduced  b>  using  the 
conformal  method  will  be  small.  On  the  other  hand,  the  added  mass  of  a  bulbous  bow  station 
may  be  significant;  thus,  when  stations  of  this  type  lie  outside  the  indicated  range,  the  close-fit 
method  should  be  used. 


2.3  Wave  Frequencies  for  Response  Computations 

Ship  response  is  computed  at  wave  frequencies  specified  in  the  program  input.  If  motions 
in  irregular  seas  are  to  be  computed,  care  must  be  taken  to  cover  the  frequency  range  of  interest 
adequately.  Physically  realistic  limits  on  the  frequency  interval  at  full  scale  are: 

0.2  <  U.V  <  2.0 


2 


where  0;^,  is  the  wave  frequency.^  SHIPMO  is  capable  of  computing  responses  for  a  maximum 
of  40  frequencies. 


2.4  Encounter  Frequencies  for  Added  Mass,  Damping,  and  Potentials 

Hull  sectional  added  mass,  damping,  and  potentials  are  calculated  at  encounter  frequencies 
specified  by  the  user.  A  maximum  number  of  40  encounter  frequencies  can  be  specified. 

The  tables  of  sectional  potentials,  added  mass,  and  damping  as  functions  of  encounter 
frequency  are  generated  (or  read  from  disk)  before  ship  motions  are  computed.  Particular  values 
are  interpolated  from  the  tables  as  necessary,  since  a  given  wave  frequency,  sea  direction,  and 
ship  speed  will  define  the  frequency  of  encounter.  Care  must  be  taken  that  the  tables  adequately 
cover  the  necessary  encounter  frequency  range.  In  general,  the  relationship  between  frequency 
of  encounter  and  the  wave  frequency  is; 

u  =  V  cos  /?,  (2.2) 

where  u  is  the  wave  encounter  frequency  and  U  is  the  ship  speed.  For  water  of  infinite  depth, 
as  assumed  in  SHIPMO,  Equation  (2.2)  can  be  expressed  as: 

u  =  Wu.  -  0.0524  F't  cos/?5  (2.3  j 

whore  Uk  is  the  ship  speed  in  knots.  It  can  be  seen,  particularly  for  a  multi-directional  analysis, 
that  a  large  encounter  frequency  range  may  exist.  In  practice,  the  added  mass,  damping,  and 
potentials  w'ill  assymptotically  approach  limiting  values  at  higher  frequencies.  A  reasonal)le  limit 
on  the  maximum  frequency  of  encounter  WEMAX  for  computations  is  WEMAX  <  20{gf 
It  is  recommended  that  the  user  check  the  output  to  ensure  that  no  significant  motions  occur 
at  encounter  frequencies  greater  than  WEM.'^X. 

In  solving  the  equations  of  motion,  a  lower  limit  WEMIN  given  as  input  is  set  on  the 
frequency  of  encounter  u;.  This  step  is  taken  because  the  linear  mathematical  model  of  ship 
motions  breaks  down  when  w  =  0  for  Wu-  >  0.  In  practice,  it  is  recommended  that  WEMIX  be 
set  such  that  WEMIN  > 


2.5  Option  to  Store  Sectional  Added  Mass,  Damping,  and  Potentials 

For  each  ship  hull  section.  SHIPMO  generates  tables  of  potentials,  added  masses,  and 
damping  over  the  range  of  encounter  frequencies  specified.  These  calculations  are  quite  time- 
consuming  and  need  not  be  repeated  in  subsequent  runs  if  the  hull  form  and  CG  location  remain 
unchanged.  Thus,  provision  is  made  that  these  tables  may  be  stored  on  disk  and  merely  read 
into  the  program  for  subsequent  runs.  This  is  most  useful  when  examining  the  positions,  types 
and  sizes  of  stabilizing  devices  for  a  single  hull. 

’Although  w  and  are  commonly  used  to  denote  wave  frequency  and  wave  encounter  frequency  respectively, 
this  manual  uses  for  wave  frequency  and  for  encounter  frequency  to  maintain  consistency  with  Reference  1. 
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3  Overview  of  Modifications 

This  section  summarizes  the  changes  made  for  version  5  of  SlIIPMO. 

3.1  Units  for  Internal  Computations 

For  previous  SHIPMO  versions,  internal  computations  are  conducted  in  FPS  units.  Interna] 
computations  for  version  5  are  done  in  metric  or  FPS  units,  depending  on  the  unit  system 
requested  for  program  output.  If  the  input  unit  system  differs  from  the  output  unit  system, 
then  input  values  are  converted  to  output  units  after  echoing  of  user  input. 

3.2  Wave  Spectra 

A  complete  list  of  wave  spectra  available  in  SHIPM05  is  given  here  for  reference. 


Quadratic  Regression  Spectrum 

The  quadratic  regression  spectrum  is  based  on  observed  wave  records  at  station  India  in 
the  North  Atlantic  [16].  Input  parameters  are  significant  wave  heiglit  11^  and  energy-averaged 
wave  period  T_i,  defined  as: 


S{T)T  AT 
J^S(T)dT 


.3.1) 


where  T  is  wave  period  and  S(T)  is  spectral  energy  density.  Tlie  spectral  energy  is 

evaluated  from  a  non-dimensional  database  as  follows: 


0  for  u-v  <  0.05  u.’_] 

*  /(//s.T^i)  for  0.05  w-_,  < 

0  for  w’u:  >  4  w'_i 


(3.2) 


where  u-'_i  —  For  fully-developed  seas,  the  peak  wave  period  Tp  used  as  input  for  the 

Bretschneider  and  JONSW.^P  spectra  can  be  related  to  the  energy  averaged  period  as  follows: 


Tp  =  1.166  T_i 


(3.3) 


Bretschneider  Spectrum 

SHIPMO  uses  the  Bretschneider  spectrum  taken  from  the  15th  Internationa]  Towing  Tank 
Conference  (ITTC)  [7]  as  follows: 


=  — — r 


-“p  “-U- 


-1948 
tT  T 

^p  '*'u  J 


(3.4) 
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JONSWAP  Spectrum 


The  JONSWAP  spectrum  models  relatively  high-peaked  spectra  typically  encountered  in 
fetch-limited  regions  [7].  The  JONSWAP  spectrum  is  obtained  by  multiplying  the  Bretschneider 
spectrum  by  a  peak  enhancement  factor  accounting  for  fetch  limited  conditions.  The  following 
JONSWAP  formulation  is  used; 


p 


,  320.^2  [-1948] 

p  L-*  p  J 

K0.159wTp-  1\^' 

—CT-)  J 


0.07  forw^  <  6.64/Tp 
0.09  for  a-v  >  6.64/Tp 


The  Bretschneider  and  JONSWAP  spectra  can  also  be  expressed  in  terms  of  the  period  Ti 
corresponding  to  the  average  frequency  using  the  following  relationship  for  fully  developed  seas: 


Ti  =  0.773  Tp 


where  Tj  is  evaluated  as  follows; 


Ti  =  2t 


SM  d  u-.. 
S(^\  )  Wu.dwu. 


Ochi  and  Hubble  Six  Parameter  Spectrum 

The  Ochi  and  Hubble  6  parameter  spectrum  models  collinear  swell  and  sea  components  as 
follows; 

4  ^  n  \  \ 

1=1  i  v^i ;  -'u- 

where  A;,  and  Wp,  are  the  spectra  shape  parameter,  significant  wave  height,  and  peak  fre¬ 
quency  for  component  i.  If  only  one  of  the  two  components  is  considered  and  the  shape  parameter 
A,  is  equals  one,  then  the  six  parameter  spectrum  is  equivalent  to  the  Bretschneider  spectrum. 


Ten  Parameter  Directional  Spectrum 

The  ten  parameter  spectrum  is  a  directional  extension  of  the  Ochi  and  Hubble  six  parameter 
spectrum.  Each  of  the  swell  and  sea  components  is  multiplied  by  its  own  directional  spreading 
function  as  follows: 


=  A(jP, )  cos' 


7  =  1. 2 
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where  u  is  the  compass  direction  (from)  of  tlie  directional  contribution,  and  P,  and  Vm,  are 
the  directional  spread  parameter  and  mean  compass  direction  (from)  for  component  i.  The 
normalization  factor  A{Pi)  is  expressed  as: 


A{P^) 


2{2P.-1)  r2(p,  +  1) 

7rr(2P.  +  l) 


i  =  1,2 


(3.12) 


Hindcast  Spectra 

Directional  spectra  from  the  Offshore  Data  Gathering  Program  (ODGP)  hindcast  model 
[10]  can  be  used  by  SHIPM05.  A  complete  description  of  hindcast  spectra  data  files  is  given  in 
Appendix  A. 4. 


User  Input  Directional  Spectra 

User  input  directional  spectra  can  be  read  by  SHIPM05  based  on  the  format  given  in 
Appendix  A.6. 


3.3  Seakeeping  Calculations 

Several  new  features  have  been  incorporated  into  the  seakeeping  calculations  at  positiojis 
on  the  ship. 


Zero-Crossing  Periods 


To  give  the  user  some  information  on  the  frequency  content  of  the  responses  in  irregular 
seas,  zero-crossing  periods,  7),,  are  now  computed  and  output.  The  zero-crossing  period  is  gi  ven 
by: 


where  m,  is  the  ith  moment  of  the  response  spectrum. 


(3.13) 


Force  Estimators  and  Motion-Induced  Interruptions 

The  new  version  of  SHIPMO  includes  an  option  for  computing  force  estimators  and  the 
resulting  incidence  of  motion-induced  interruptions  (Mil’s)  [5].  The  lateral  and  longitudinal  force 
estimators  indicate  the  force  due  to  ship  motion  acting  on  a  person  or  object  relative  to  the  local 
ship  axes.  Retaining  linear  terms  and  expressing  as  force  per  unit  mass,  the  longitudinal,  lateral, 
and  vertical  force  estimators  are  as  follows: 


Plong  — 

-0, 

+  9  ^5 

(3.14) 

Plat  - 

-02 

-  9  'h 

(3.1.5) 

Pv  = 

-«3 

-  9 

(3.1C) 
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where  oi,  02  and  03  are  the  local  accelerations  in  the  x,  y  and  z  directions,  and  1J4  and  7/5  are 
the  ship  roll  and  pitch  angles  respectively. 

In  order  to  determine  the  effects  of  the  above  forces  on  equipment  or  people  aboard  ship, 
Graham  [5]  introduced  a  frequency-domain  method  for  estimating  the  incidence  of  equipment 
sliding,  or  of  personnel  loss-of- balance  events,  termed  motion-induced  interruptions. 

First  consider  an  object  on  deck  and  suppose  that  the  object  has  a  static  coefficient  of 
friction  given  by  /r,.  In  this  case,  slides  can  only  occur  in  the  port  or  starboard  directions.  A 
slide  to  port  will  occur  whenever: 


^lat  ^  F y 


(3.17) 


or 


-  02  -  5  r?4  -  //j  03  > 


(3.18) 


The  quantity  on  the  left  hand  sides  of  Equation  (3.18)  was  called  a  generalized  lateral  force 
estimator  in  Reference  5.  In  the  present  work,  this  quantity  will  be  called  the  port  sliding 
estimator  function,  and  denoted  by  Sport-  Tipping  estimator  functions  will  be  considered  below. 

In  Reference  5,  the  number  of  slides  per  unit  time  was  derived,  under  the  assumption  that 
Sport{t)  followed  the  Rayleigh  distribution.  In  fact,  the  result  holds  provided  that  SportlO  is 
a  zero- mean  Gaussian  process,  and  is  not  restricted  to  narrow-band  processes.  Let 
denote  the  (one-sided)  power  spectral  density  of  SportiO  ^  ^  function  of  encounter  frequency  w, 
and  denote  by  m„  the  spectral  moments  of  55p^„(u;): 


w„ 


(3.19) 


Rice  [17]  (see  also  [18])  showed  that  the  number  of  upcrossings  (or  downcrossings)  per  unit  time 
of  a  threshold  at  level  fi,g  is  given  by: 


^iport 


1  (//5  gf 

2  TT  Y  mo  2  mo 


(3.20) 


The  quantity  Mport  is  the  number  of  slides  per  second  in  the  pert  direction.  For  a  shipboard 
task  requiring  Tj  seconds  to  complete,  the  expected  number  of  port  sliding  incidents  during  the 
duration  of  the  task  is  obtained  by  multiplying  Equation  (3.20)  by  Tj- 

Similarly,  the  number  of  sliding  incidents  that  occur  in  the  starboard  direction  can  be 
determined  from  the  number  of  upcrossings  of  the  following  sliding  estimator  function: 


Ssibd  —  ^^2  +  J74  ~  A's  03 


(3.21) 


Now  consider  the  problem  of  motion-induced  interruptions.  First  consider  a  person  facing 
in  the  fore-and-aft  direction,  and  assume  that  his  centre  of  gravity  is  located  at  a  height  h  above 
the  deck,  and  that  his  stance  width  is  2s,  as  shown  in  Figure  5.  By  taking  moments  about  the 
lef'  oot  position,  it  can  be  seen  that  the  person  will  tip  to  port  if: 

-02-994-^03  >  ^g  (3.22) 


i 


The  quantity  on  the  left  hand  side  of  Equation  (3.22)  will  be  called  the  port  tipping  estimator 
function,  and  will  be  denoted  by  Tport-  As  in  Reference  5,  the  quantity  s/h  is  referred  to  as 
the  tipping  coefficient.  The  number  of  port  tipping  incidents  can  be  determined  as  in  Equa¬ 
tion  (3.20).  From  typical  personnel  stance  geometries,  representative  values  for  s/h  would  be 
0.25  and  0.17  for  a  person  facing  forward  and  sideways,  respectively.  On  the  other  hand,  a 
typical  static  coefficient  for  dry  deck  conditions  is  0.7.  Comparing  Equations  3.18  and  3.22,  it  is 
clear  that  for  dry  deck  conditions  and  non-skid  surfaces,  people  will  tip  over  before  they  slide. 

Now  consider  tipping  in  the  presence  of  combined  vertical,  lateral,  and  longitudinal  ac¬ 
celerations.  Suppose  that  a  person  is  facing  forward,  and  that  at  a  given  instant  in  time  the 
resultant  force  parallel  to  the  deck  (the  vector  sum  of  Fiat  and  Fiong)  is  in  the  direction  6  shown 
in  Figure  6.  The  corners  of  the  rectangle  shown  in  this  figure  are  defined  by  the  location  of  the 
feet  of  the  crew  member.  The  resultant  force  will  result  in  a  tip  in  the  direction  6  provided  that 


nong  COS  6  +  Flat  Sin  d  -  -j- 


(3.23) 


For  a  tip  to  occur  at  some  point  Os  on  the  top  of  the  rectangle,  it  is  necessary  and  sufficient 
that 


Flo 


m 


0.3  > 


/(O) 


long  ^  ^ 

To  see  this,  note  that  Equation  3.23  is  equivalent  to 


/fL,  +  n..  -  X  ^ 


/(^) 


Multiplying  this  result  by  cos^  =  Fiong! sj Ff^^^  -|-  we  obtain 


(3.24) 


(3.25) 


l{6)  cosf  l{6)  cos 6 

Fiong  -  ^  03  >  9 

which  is  equivalent  to  Equation  3.24. 

If  the  forward  tipping  estimator  function  is  defined  as 


^  fore 


-a\  +  9  m  -  —  03 


(3.26) 


(3.27) 


then  the  forward  tips  can  be  counted  by  determining  the  uperossings  of  7/ore-  Similarly,  the  aft 
tipping  estimator  function  can  be  defined  via 


Tajt  =  ai  -  grj5 - —03  (3.28) 

The  total  number  of  motion-induced  interruptions  is  determined  by  summing  the  uperossings 
of  the  four  tipping  estimator  functions. 

SHIPM05  requires  the  user  to  input  the  duration  of  the  task  of  interest,  T7-.  plus  values 
for  the  lateral  and  longitudinal  tipping  coefficients.  Representative  values  for  these  coeffiicients 
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Table  1;  Representative  Tipping  Coefficients 


Tipping  forward  0.17 

Tipping  sideways  0.25 

Sliding  0.7 


Table  2:  Mil  Risk  Levels 
Risk  Level  Mil’s  per  Minute 


1. 

possible 

0.1 

2. 

probable 

0.5 

3. 

serious 

1.5 

4. 

severe 

3.0 

5. 

extreme 

5.0 

given  above  are  also  given  again  in  Table  1.  As  a  guide  to  the  user  in  interpreting  the  incidence 
of  Mil  events,  Table  2  reproduces  information  from  Reference  5. 

In  the  case  of  negligible  longitudinal  acceleration,  which  is  usual  on  frigates  and  destroyers, 
the  user  can  determine  the  sliding  incidence  by  inputting  the  static  coefficient  of  friction  for 
both  the  lateral  and  longitudinal  tipping  coefficients.  The  user  is  cautioned  that  in  the  presence 
of  significant  longitudinal  acceleration  it  is  not  true  that  the  total  number  of  slides  is  the  sum 
of  the  uperossings  of  the  sliding  estimator  functions  in  the  forward,  port,  aft.  and  starboard 
directions.  The  reason  for  this  is  that  the  coefficient  of  friction  is  the  same  in  all  directions, 
which  implies  that  slides  can  occur  in  oblique  directions  without  an  associated  uperossing  by 
any  of  the  sliding  estimator  functions  iii  the  four  primary  directions. 


Off-Centre  Seakeeping  Positions 


Version  5  of  SHIPMO  reads  y  coordinates  for  seakeeping  positions  in  addition  to  x  and  x 
coordinates.  If  a  seakeeping  position  is  located  off  the  ship  centreline,  then  responses  for  headings 
180-360  degrees  are  not  symmetrical  to  responses  for  headings  0-180  degrees.  Consequently, 
input  sea  directions  should  range  from  0  to  at  least  315  degrees  if  any  seakeeping  position  is 
located  off-centre. 
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Propeller  Emergence  Calculations 

SHIPM05  includes  an  option  for  predicting  propeller  motions  and  emergence,  which  can 
influence  added  power  in  waves.  The  theory  described  in  Reference  19  for  predicting  keel  emer¬ 
gence  is  used  for  predicting  propeller  emergence.  The  propeller  emergence  calculations  can 
include  the  effects  of  the  dynamic  waterline  due  to  ship  forward  motion  and  wave  diffraction 
by  the  hull;  however,  the  empirical  methods  for  predicting  the  dynamic  waterline  may  give 
inaccurate  results  at  the  stern. 

3.4  Hydrodynamic  End-Effect  Terms 

The  new  version  of  SHIPMO  has  an  option  for  computing  hydrodynamic  end-effect  terms 
[13]  for  transom  stern  ships  at  forward  speed.  End-effect  terms  are  expected  to  give  somewhat 
improved  motion  predictions  at  higher  forward  speed  {Fn  >  0.2)  [14]. 

3.5  Surge  Excitation  Forces 

A  new  method  has  been  implemented  for  computing  surge  excitation  force.  An  average 
surge  normal  for  a  finite  length  ship  section  is  determined  by  considering  the  hull  area  projected 
onto  the  y  —  z  plane.  Adopting  the  notation  of  Reference  1  and  assuming  that  wavelength  is 
large  relative  to  draft  and  beam,  the  sectional  and  total  surge  excitation  force  are; 


/i(x)  =  5  exp(-t  Atu.  2:  COS/3,) 

dx 

(3.29) 

Fi  =  pa  j^fi{x)dx 

(3.30) 

where  Ai  is  sectional  area,  p  is  water  density,  and  a  is  wave  amplitude.  The  contribution  from 
the  transom  to  the  total  excitation  force  is  done  by  taking  the  limit  dx  —  0  in  Equations  3.29 
and  3.30.  The  transom  contribution  to  surge  excitation  is  included  regardless  of  whether  the 
end-effect  terms  of  the  previous  section  are  requested. 


4  Input  Description 

This  section  describes  the  main  features  of  input  for  SHIPM05.  Detailed  record  descrip¬ 
tions  and  a  sample  input  file  are  given  in  Appendix  A. 

4.1  General  Input  Format 

Program  input  consists  of  an  alphanumeric  title  and  records  of  free-format  numerical  data. 

Two  systems  of  units  are  available  in  SHIPM05,  either  metric  or  British,  as  shown  in 
Table  3.  Most  internal  computations  are  carried  out  in  British  units;  however,  some  R.M.S 
computations  are  carried  out  in  the  .system  of  units  selected  for  output. 
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Table  3:  Units  for  Metric  and  British  Systems 


Metric 

British 

Length 

metres 

feet 

Displacement 

tonnes 

long  tons 

Force 

Newtons 

pounds 

Speed 

knots 

knots 

A  sample  input  file  is  given  in  Appendix  A.3  for  a  frigate.  The  principal  particulars  of  the 
vessel  are  listed  in  Table  4.  A  body  plan  is  shown  in  Figure  7.  The  vessel  has  twin  screws,  twin 
rudders,  and  a  pair  of  bilge  keels. 


Table  4:  Ship  Particulars  for  Sample  Input  File 


Displacement,  A 

2932  tonnes 

Length  between  perpendiculars,  L 

108.5  m 

Block  coefficient,  Cb 

0.50 

Length/beam,  LjB 

8.52 

Length/draft,  LfT 

2C.0 

Metacentric  height /beam.  CM  IB 

0.05 

Roll  radius  of  gyration/beam  (in  water),  kx^lB 

0.35 

4.2  Changes  to  Input  from  Previous  SHIPMO  Versions 

When  creating  new  input  files  for  Version  5  of  SHIPMO,  attention  should  be  given  to 
changes  in  input  due  to  new  program  features.  The  reader  should  also  note  that  most  vertical 
coordinates  are  now  referenced  to  the  baseline  to  make  input  files  more  consistent.  A  complete 
list  of  input  changes  is  given  in  Appendix  A. 2. 


4.3  Running  Input  Files  for  Previous  Versions  of  SHIPMO 

Although  the  above  subsection  indicates  that  input  for  SHIPM05  differs  from  previous 
versions  of  SHIPMO,  version  5  is  capable  of  reading  input  files  for  old  versions  without  modifi¬ 
cation.  If  the  new  variable  lEND  is  not  present  in  Record  (d)  of  an  input  file,  then  the  input  file 
is  assumed  to  be  for  earlier  versions  of  SHIPMO.  Modification  of  input  variables  is  subsequently 
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done  within  SHIPM05.  The  only  possible  source  of  incompatibility  in  old  input  files  is  the  sea 
direction  PSDIR(I)  for  short-crested  sea  computations,  for  which  a  consistency  check  has  been 
added  to  SHIPM05. 


5  Output  Description 

This  section  gives  a  brief  summary  of  program  output.  A  complete  description  of  program 
output  and  sample  output  is  given  in  Appendix  B. 

5.1  Printer  Output 

Main  output  for  SHIPMO  is  written  to  an  ASCII  file  suitable  for  printing  on  a  line  printer. 
A  description  of  program  output  is  given  in  Appendix  B.l.  A  sample  output  file  is  given  in 
Appendix  B.3. 

5.2  Post-Processor  File 

If  post-processing  output  is  requested  from  SHIPMOo  for  use  in  other  programs,  an  un¬ 
formatted  file  is  written  to  disk.  Post-processing  output  has  been  updated  from  the  previous 
output  described  in  Reference  1,  and  is  given  in  Appendix  B.2. 

6  Concluding  Remarks 

This  user's  manual  has  described  the  improvements  implemented  in  the  latest  version  of 
ship  motion  program  SHIPMO,  which  has  been  designated  SHIPM05.  These  improvements  are 
primarily  in  three  areas.  First,  an  enhanced  capability  for  wave  modelling  has  been  incorporated 
in  the  projram  with  the  addition  of  JONSWAP,  10-parameter,  and  hindcast  and  user  input 
directional  wave  spectra.  Second,  the  program  now  computes  forces  relative  to  local  ship  axes 
and  the  resulting  incidence  of  motion-induced  interruptions  or  sliding  events  at  user-specified 
seakeeping  positions,  which  are  no  longer  restricted  to  the  centreline  of  the  ship.  Finally,  motion 
computations  may  now  include  hydrodynamic  end-effects. 

All  input  vertical  coordinates  are  referenced  to  the  ship  baseline  in  the  new  version  of  the 
program.  Post-processor  output  has  also  been  updated. 

Example  input  and  output  files  are  presented  for  the  updated  program.  Details  of  the  input 
and  output  record  structures  are  presented  in  the  appendices. 


12 


T^i  =  SURGE  7j3=  HEAVE  7^5=  PITCH 

7^2=  SWAY  7J4=R0LL  7^6=  yaw 

Figure  1;  Axis  System  (Arrows  Indicate  Senses  of  Motion) 


y 


Figure  2:  Definition  of  Sea  Direction  Relative  to  Ship  Heading 
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Figure  3:  Definition  of  Sea  Direction  and  Ship  Heading  Relative  to  Compos 


0.5  1.0  1.5  2.0 


HALF  BEAM,  DRAFT  RATIO.  Bx  /2Tx 


Figure  4:  Range  of  Sectional  Area  Coefficients  for  Lewis  and  MIT  Forms 


IT 


Figure  6:  Tipping  Geometry  in  the  Presence  of  Combined  Vertical,  Lateral,  and  Longitudinal 
Forces 


LOAD  WATERLINE 


A- BRACKET 


Figure  13:  Foil  Inputs 


Figure  14:  Tank  Inputs 


Appendix 

A  Program  Input 

A.l  Input  Records 

Detailed  descriptions  of  SHIPM05  input  records  are  given  below.  Each  new  input  record 
or  sub-record  corresponds  to  a  new  file  line.  The  format  of  the  input  file  may  be  adjusted 
by  inserting  extra  blanks  between  any  numerical  data,  and  by  placing  data  from  within  any 
particular  record  on  separate  lines;  however,  separate  records  cannot  be  combined  on  a  single 
line. 

Record  (a),  Fifty  Character  Title 
TITLE  (columns  1  -  80) 

TITLE  Alphanumeric  title  (maximum  of  80  characters)  which  is  written  on  output. 

Record  (b),  Control  Integers 

IN,  lOUT,  IRHO,  ISPEC,  ICORR,  IRESP,  IPPF  (7  integers) 

IN  Specifies  the  system  of  units  used  for  input  data. 

IN  =  0  British  units. 

IN  =  1  Metric  units. 

lOUT  Specifies  the  system  of  units  used  for  output  data. 

lOUT  =  0  Britisli  units. 
lOUT  =  1  Metric  units. 

IRHO  Specifies  whether  ship  is  in  salt  or  fresh  water. 

IRHO  =  0  Salt  water  (p  =  1.9906  slug/ft^,  1026.08  kg/m^). 

IRHO  =  1  Fresh  water  (p  =  1.9-40  slug/ft^,  1000.0  kg/m^). 

ISPEC  Specifies  the  seaway  spectrum  to  be  used  for  motion  calculations  in  irregular 
seaways. 

ISPEC  =  0  Quadratic  regression  spectrum. 

ISPEC  =  1  Bret.schneider  two  parameter  spectrum. 

ISPEC  =  2  Measured  uni-directional  spectrum. 

ISPEC  =  3  JONS\V.4P  spectrum. 

ISPEC  =  4  Ochi  and  Hubble  six  parameter  uni-directional  spectrum. 
ISPEC  =  5  Ten  parameter  directional  spectrum. 
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ISPEC  =  6  ODGP  hindcast  directional  spectrum. 

ISPEC  =  7  User  input  directional  spectrum. 

ICORR  Specifies  corrections  to  be  made  to  the  incoming  wave  train. 

ICORR  =  0  No  corrections  are  applied. 

ICORR  =  1  Corrections  are  made  for  wave  deformation  by  the  hull  (dy¬ 
namic  swell-up)  and  for  wave  profile  generated  from  hull  for¬ 
ward  motion. 

ICORR  =  2  A  correction  is  made  only  for  wave  profile  generated  from  hull 
forward  motion. 

IRESP  Control  integer  for  program  output  of  regular  wave  responses. 

IRESP  =  0  Suppress  output  of  regular  wave  responses,  and  of  roll  damping 
coefficients. 

IRESP  =  1  Allow  output  of  regular  wave  responses,  but  suppress  output 
of  roll  damping  coefficients. 

IRESP  =  2  Allow  output  of  regular  wave  responses  and  of  roll  damping 
coefficients. 

IPPF  Control  integer  governing  disk  file  storage  of  resulting  frequency  responses  and 
RMS  motions  for  post-processing. 

IPPF  =  0  Disk  storage  is  not  used. 

IPPF  =  1  Frequency  responses  and  RMS  motions  are  stored  in  a  sequen¬ 
tial  disk  file  for  post-processing.  The  name  of  the  disk  file  must 
be  specified  in  Record  (bl).  File  format  is  given  in  Appendi.x 
B.2. 

Note:  If  a  measured  seaway  uni-directional  spectrum  is  being  used  (ISPEC  =  2)  the 
units  of  the  spectral  density  must  be  specified  in  Record  (f2)  below.  The  units 
of  the  spectral  density  may  be  chosen  independently  of  the  value  of  IN  specified. 
In  the  case  of  a  hindcast  spectrum  (ISPEC  =  6)  the  input  is  always  in  metric 
units,  as  described  in  Appendix  A.4. 


Record  (bl),  File  Name  if  Disk  Storage  Specified  for  Responses  and  Motions 
Read  only  if  IPPF  =  1  in  Record  (b) 

PPNAME  (columns  1  -  20) 

PPNAME  The  disk  file  name  for  storing  frequency  responses  and  RMS  motions  for  post¬ 
processing  (maximum  20  characters). 
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Record  (c),  Wave  Frequencies  for  Response  Calculations 
WMIN,  WMAX,  D\V  (3  reals) 

See  discussion  in  Section  2.3 

WMIN  Lowest  wave  frequency  for  which  ship  motions  are  to  be  calculated  (rad/sec). 

WMAX  Highest  wave  frequency  for  which  ship  motions  are  to  be  calculated  (rad/sec). 

DW  Increment  in  wave  frequency  between  WMIN  and  WMAX  for  which  ship  motions 
are  to  be  calculated  (rad/sec).  Due  to  storage  considerations,  DW  should  be  set 
such  that  DW  >  (WMAX  -  WMlN)/39. 

Record  (d),  Encounter  FVequencies  for  Added  Mass  and  Damping  Calculations 
ISAVE,  WEMIN,  WEMAX,  DWE,  METHOD,  lEND  (1  integer,  3  reals,  2  integers) 

See  discussion  in  Sections  2.4  and  2.5. 

ISAVE  Control  integer  governing  disk  file  storage  of  the  two-dimensional  section  poten¬ 
tials,  added  mass  and  damping. 

ISAVE  =  0  Disk  storage  is  not  used. 

ISAVE  =  1  The  hull  sectional  potentials,  added  mass  and  damping  which 
are  computed  in  the  program  are  stored  in  a  sequential  disk 
file  for  future  use.  The  name  of  the  disk  file  must  be  specified 
in  Record  (dl). 

ISAVE  >  2  The  hull  sectional  potentials,  added  mass  and  damping  for  the 
ship  are  read  from  a  sequential  disk  file  which  has  already  been 
created  through  a  previous  execution  of  the  program.  The  disk 
file  name  must  be  specified  in  Record  (dl). 

WEMIN  Lowest  encounter  frequency  for  which  the  hull  sectional  potentials,  added  mass 
and  damping  are  calculated  (rad/sec). 

WEMAX  Highest  encounter  frequency  for  which  the  hull  sectional  potentials,  added  mass 
and  damping  are  calculated  (rad/soc). 

DWE  Increment  in  encounter  frequency  between  WEMIN  and  WEMAX  for  which  the 
hull  sectional  properties  are  calculated  (rad/sec).  Due  to  storage  considerations, 
DWE  should  be  set  such  that  DWE  >  (WEMAX  -  WEMIN)/39. 

METHOD  Control  integer  governing  the  computational  method  for  the  hull  sectional  prop¬ 
erties.  See  discussion  in  Section  2.2. 

METHOD  =  0  Hull  sectional  properties  are  computed  using  section  offsets 
and  the  close-fit  method.  No  long  wavelength  approximations 
are  made  in  computing  exciting  forces. 

METHOD  =  1  The  close-fit  method  is  used  as  for  METHOD  =  0;  however, 
long  wavelength  approximations  are  made  in  computing  excit¬ 
ing  forces. 


METHOD  =  2  Hull  sectional  properties  are  computed  using  general  section 
data  (beam,  draft,  area  coefficient)  and  the  conformal  map¬ 
ping  method.  Long  wavelength  approximations  are  made  in 
computing  exciting  forces.  This  method  involves  less  compu¬ 
tation  time. 

lEND  Control  integer  for  end-effect  hydrodynamic  terms. 

lEND  =  0  Hydrodynamic  end-effect  terms  not  included. 
lEND  =  1  Hydrodynamic  end-effect  terms  included. 

Record  (dl),  File  Name  for  Storing  Potentials,  Added  Mass  and  Damping 
Read  only  if  ISAVE  >  0  in  Record  (d). 

ADNAME  (columns  1  -  20) 

ADNAME  The  disk  file  name  for  hull  sectional  potentials,  added  mass  and  damping  coeffi¬ 
cients  (maximum  20  characters). 

Record  (e),  Number  of  Sea  Directions  and  Short-Crested  Analysis  Spreading  Angle 
Read  only  if  ISPEC  <  4  in  Record  (b). 

NSD,  ANGLE  (1  integer,  1  real) 

NSD  Number  of  principal  sea  directions  to  be  considered  with  respect  to  ship  heading 
(maximum  of  36).  If  NSD  =  0,  default  headings  from  0  to  180  degrees  in  15 
degree  increments  and  a  spreading  angle  of  90  degrees  are  used. 

ANGLE  The  spreading  angle  (degrees)  to  be  used  in  a  short-crested  sea  spectrum  analysis. 
If  ANGLE  <  0.0  and  NSD  >  0,  no  short-crested  analysis  is  carried  out. 

Record  (el),  Principal  Sea  Directions 

Read  only  if  ISPEC  <  4  in  Record  (b)  and  NSD  >  0  in  Record  (e). 

PSDIR(I)  (NSD  reals) 

PSDIR(I)  The  principal  sea  directions  to  be  considered  relative  to  the  ship  velocity  vector 
(degrees).  There  must  be  NSD  values. 

Notes:  1.  If  ANGLE  >  0.0  and  NSD  >  0  in  Record  (e),  a  sufficient  range  of  sea 

directions  must  be  given  for  short-crested  calculations.  For  all  seakeeping 
positions  on  the  centreline  (YPOS(I)  =  0.0  for  all  I  in  Record  (jl)),  PSDIR 
must  cover  the  range  0  -  180  degrees.  If  any  seakeeping  position  is  off 
centreline,  PSDIR  must  cover  the  range  from  0  to  a  minimum  of  315.0 
degrees. 

2.  If  ANGLE  >  0.0  and  NSD  >  0  in  Record  (e),  PSDIR  values  must  be  given  in 
ascending  or  descending  order,  with  a  maximum  increment  of  45.0  degrees 
between  any  2  values. 


24 


Record  (f),  Sea  State  Definition 
NSEA  (1  integer) 

NSEA  Number  of  seaways  for  which  motions  are  to  be  computed  (maximum  10). 

Notes:  1.  If  motions  in  irregular  seas  are  not  desired,  set  NSEA  =  0.  Record  (fl)  is 

required. 

2.  If  a  measured  (ISPEC  =  2),  hindcast  (ISPEC  =  6),  or  user  input  directional 
(ISPEC  =  7)  spectrum  is  used,  NSEA  must  be  set  equal  to  1. 


Record  (fl),  Wave  Slope  if  no  Sea  States  Specified 
Read  only  if  NSEA  =  0  in  Record  (f). 

WVSLP  (1  real) 

WVSLP  A  wave  slope  (in  degrees)  used  to  determine  wave  height  in  the  estimation  of 
viscous  roll  damping.  A  nominal  slope  of  3  degrees  is  recommended. 


Record  (f2),  Significant  Wave  Height  and  Characteristic  Period 

Read  only  if  ISPEC  <  3  in  Record  (b)  and  NSE.4  >  0  in  Record  (f). 
HSW(I),  TSW(I)  (NSEA  records,  2  reals/record) 


HSW(I)  Significant  wave  height  of  seaway  I  (ft  or  m). 
TS\V(I)  Characteristic  wave  period  (sec): 


ISPEC  =  0 
ISPEC  =  1 
ISPEC  =  2 
ISPEC  =  3 


energy- averaged  period  of  seaway  I. 

peak  wave  period  of  seaway  I. 

any  measure  of  the  wave  period  may  be  used. 

peak  wave  period  of  seaway  I. 


Note:  For  ISPEC  =  2,  the  significant  wave  height  and  period  are  used  to  label  tables  in 
the  output,  and  HSW  is  also  used  to  calculate  the  sea  state  number.  HSW  and 
TSW  do  not  affect  the  response  calculations.  If  HSW  and  TSW  are  not  known 
exactly,  estimates  may  be  used. 


Record  (f3),  Input  of  Measured  Point  Spectrum 

Read  only  if  ISPEC  =  2  in  Record  (b)  and  NSEA  =  1  in  Record  (f). 
ITOP,  ISEAIN,  INDATA  (3  integers,  free  format) 

ITOP  Number  of  frequencies  for  spectral  density  input  (maximum  300). 
ISEAIN  System  of  units  for  spectral  density: 


ISEAIN  =  0  British  units 


ISEAIN  =  1  Metric  units 


The  units  of  spectral  density  (m^/(rad/s)  or  ft^/(rad/s))  may  be  chosen  inde¬ 
pendently  of  the  value  IN  specified  in  Record  (b). 

IN  DATA  Format  of  the  spectral  density  input 

INDATA  =  0  Spectral  data  in  the  form  of  output  from  the  DREA  WAVE 
program.  SHIPMO  calls  a  file  named  WAVE.LST  for  the 
data,  which  is  echoed  in  the  SHIPM05  output.  A  typical 
WAVE.LST  file  is  shown  in  Appendix  B  of  the  SHIPM04  man¬ 
ual  [4]. 

INDATA  =  1  Spectral  data  is  stored  in  a  file,  the  name  of  which  must  be 
specified  in  Record  (f3a). 

Record  (fSa),  File  Name  for  Spectral  Data  from  Disk  File 
Read  only  if  INDATA  =  1  in  Record  (f3). 

SPECTR  (columns  1  -  20) 

SPECTR  The  disk  file  name  for  wave  sea  spectral  data  (maximum  20  characters).  The 
spectral  data  file  must  contain  2  x  ITOP  entries  in  the  following  order; 

^u'l  1  •S'('*'’u.l  )i  •S^('*'V'2)*^’u;3?  ■S^(“'V-3)i  ^wITOPi  S{u^uiITOP) 

where  is  the  spectral  density  at  frequency  uv,.  A  sample  input  file  is 

given  in  Appendix  C  of  the  SHIPM04  manual  [4], 


Record  (f4),  Ochi  and  Hubble  6  Parameter  Uni-Directional  Spectrum 
Read  only  if  ISPEC  =  4  in  Record  (b)  and  NSE.4  >  0  in  Record  (f). 

ZETAl(I),  OMEGAl(I).  LAMBl(I),  ZET.42(I).  OMEGA2(I).  LAMB2(1)  (NSEA 
records,  6  reals/record) 


ZETAl(I) 
OMEGA  1(1) 
LAMBl(I) 
ZETA2(I) 
OMEGA2(I) 
LAMB2(I) 


Significant  wave  height  of  swell  component  (ft  or  m). 
Peak  frequency  of  swell  component  (rad/s). 

Spectral  shape  parameter  of  swell  component. 
Significant  wave  height  of  sea  component  (ft  or  m). 
Peak  frequency  of  sea  component  (rad/s). 

Spectral  shape  parameter  of  sea  component. 
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Record  (f5),  10  Parameter  Directional  Spectrum 

Read  only  if  ISPEC  =  5  in  Record  (b)  and  NSEA  >  0  in  Record  (f). 

ZETAl(I),  OMEGAl(I),  LAMBl(I),  KAPAl(I),  EXPPl(I),  ZETA2(I), 
0MEGA2(1),  LAMB2(1),  KAPA2(I),  EXPP2(I)  (NSEA  records,  10  reals/record) 

ZETAl(I)  Significant  wave  height  of  swell  component  (ft  or  m). 

OMEGAl(I)  Peak  frequency  of  swell  component  (rad/s). 

LAMB  1(1)  Spectral  shape  parameter  of  swell  component. 

KAPAl(I)  Mean  compass  direction  (from)  of  swell  component  (degrees). 

EXPPl(I)  Spreading  parameter  of  swell  component. 

ZETA2(I)  Significant  wave  height  of  sea  component  (ft  or  m). 

OMEGA2(l)  Peak  frequency  of  sea  component  (rad/s). 

LAMB2(I)  Spectral  shape  parameter  of  sea  component. 

KAPA2(1)  Mean  compass  direction  (from)  of  sea  component  (degrees). 

EXPP2(I)  Spreading  parameter  of  sea  component. 


Record  (f8),  Hindcast  Spectrum  File  Name 

Read  only  if  ISPEC  =  6  in  Record  (b)  and  NSEA  =  1  in  Record  (f). 

HCFILE  (columns  1  -  20) 

HCFILE  The  disk  file  name  from  which  hindcast  data  are  read  (maximum  20  characters). 

The  fo"mat  of  hindcast  spectral  data  files  is  described  in  Appendix  A. 4.  The 
units  for  this  input  file  are  independent  of  the  values  chosen  for  IN  and  lOUT  in 
Record  (b). 


Record  (f7),  User  Input  Directional  Spectrum  File  Name 

Read  only  if  ISPEC  =  7  in  Record  (b)  and  NSE.4  =  1  in  Record  (f). 

SPDIRFILE  (columns  1  -  20) 

SPDIRFILE  The  disk  file  name  from  which  the  user  input  directional  spectrum  is  read  (max¬ 
imum  20  characters).  The  format  of  the  user  input  directional  spectrum  files  is 
described  in  Appendix  A.C.  The  units  for  this  input  file  are  independent  of  the 
values  chosen  for  IN  and  lOUT  in  Record  (b). 


Record  (g),  Ship  Speeds 

UKMIN,  UKMAX,  DUK  (3  reals) 

UKMIN  Lowest  ship  speed  for  which  motions  are  to  be  calculated  (knots). 
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UKMAX  Highest  ship  speed  for  wliich  motions  arc  to  be  calculated  (knots). 

DUK  Increment  in  ship  speed  between  I’KMIN  and  UKMAX  for  which  ship  motions 
are  to  be  calculated  (knots). 

Record  (h),  Basic  Ship  Data 

EL,  KG,  GMIN,  ZWL,  RNF,  RRGB,  YRGL,  IRG  (7  reals,  1  integer) 

EL  Length  between  perpendiculars  (ft  or  m). 

KG  Height  of  CG  (ft  or  m)  above  the  baseline. 

GMIN  Metacentric  height  (ft  or  m).  If  input  value  is  0.0,  GMIN  will  be  computed  from 
offset  data  and  will  not  be  corrected  for  internal  free  surfaces. 

ZWL  Vertical  position  (z  coordinate,  ft  or  m)  of  the  waterplane  with  respect  to  the 
datum  line  used  for  hull  and  appendage  offsets.  Only  used  if  offsets  are  being 
input  up  to  the  load  waterline  (i.e.  NOFF  <  8  for  all  stations  in  Record  (il) 
below).  Otherwise,  input  ZWL  =  0.0. 

RNF  Roll  natural  frequency  (rad/sec).  Used  if  RRGB  =  0.0. 

RRGB  Roll  radius  of  gy  ration,  non-dimensionalized  with  respect  to  beam.  Used  if  RNF 

=  0.0. 

YRGL  Yaw  radius  of  gyration,  non-dimensionalized  with  respect  to  length  (EL).  Also 
used  for  pitch  radius  of  gyration.  If  unknown.  0.25  can  be  used  as  a  representative 
value. 

IRG  Control  integer  for  roll  radius  of  gyration  calculation: 

IRG  =  0  1-1  =  A(RRGB  x  B)^  -  RRGB  or  RNF  has  been  measured 
in  water.  RRGB  lies  between  0.35  and  0.40  for  most  small 
warships. 

IRG  =  1  14  =  A(RRGB  x  B)^ .  RRGB  or  RNF  has  been  measured  in 

air. 

Note;  One  of  RNF  or  RRGB  must  be  specified,  with  the  other  input  as  0.0. 

Record  (i),  Station  Definition  and  Scaling  Factors 

NSTOT,  NST,  BMSF,  DTSF  (2  integers,  2  reals) 

NSTOT  Total  number  of  equally  spaced  stations  into  which  the  ship  hull  has  been  divided 
for  representation  (ma.ximum  21). 

NST  The  number  of  stations  for  which  offsets  or  beam,  draft  and  area  coefficient  are 
input.  NST  can  be  le.ss  than  .NSTOT  if  the  ship  is  lacking  a  bulbous  bow  or 
transom  stern. 


BMSF  A  beam-wise  scaling  factor,  applied  to  the  horizontal  offsets,  YA,  given  in  Record 
(i2a),  or  to  the  station  beam  given  in  Record  (i3).  A  default  value  of  1.0  is 
assigned  if  BMSF  <  0. 

DTSF  A  draft-wise  scaling  factor,  applied  to  the  vertical  offsets,  ZA,  given  in  Record 
(i2b),  or  to  the  station  draft  given  in  Record  (i3).  A  default  value  of  1.0  is 
assigned  if  DTSF  <  0. 

Note:  NST  should  be  sufficiently  large  to  adequately  describe  the  hull.  For  each  of 
the  NST  stations,  there  must  be  one  Record  (il)  followed  by  either  one  each  of 
Records  (i2a)  and  (i2b),  or  one  Record  (i3). 

Record  (il),  Station  Number,  Number  of  Offsets  and  Eddy-Making  Calculations 
XA(1),  NOFF(I)  lEDDY(l)  (1  real,  2  integers) 

XA(I)  Station  number  (0  <  XA  <  NSTOT  -  1).  Station  0  is  at  the  forward  perpen¬ 
dicular,  and  station  NSTOT  —  1  is  at  the  after  perpendicular. 

NOFF(l)  Number  of  offsets  to  follow  in  Records  (i2a)  and  (i2b)  (maximum  12). 

NOFF(I)  <  1  Beam,  draft,  and  area  coefficient  will  be  input  (Record  (i3))  to 
describe  the  hull  section,  which  will  be  modelled  using  a  Lewis 
or  M.l.T.  form.  The  range  of  section  parameters  for  which 
Lewis  and  M.l.T.  forms  exist  is  shown  in  Figure  4  which  has 
been  adapted  from  Reference  15. 

NOFF(I)  >  1  Offsets  will  be  input  (Records  (i2a)  and  (i2b)  up  to  the  load 
waterline  or  deck  edge.  If  NOFF(I)  <  8  for  all  I,  offsets  are 
input  up  to  the  load  waterline.  If  NOFF(I)  >  8  for  some  1,  then 
offsets  are  input  up  to  the  deck  edge  for  all  stations,  and  the 
program  will  calculate  the  load  waterline  location.  In  this  case. 
Records  (i4),  (i5)  and  (i6)  below  specify  the  loading  conditions. 

lEDDY(I)  Control  integer  for  eddy-making  calculations. 

lEDDY  =  0  The  station  has  a  shape  unlikely  to  produce  eddies  as  the  ship 
rolls,  such  as  a  destroyer  hull  section  with  extremely  rounded 
bilges. 

lEDDY  =  1  For  V  or  U-shaped  stations  such  as  normally  occur  far  forward. 

lEDDY  =  2  For  stations  triangular  at  the  keel,  or  stations  spanned  by  the 
skcg, 

lEDDY  =  3  For  stations  with  a  shape  likely  to  produce  eddies  at  the  bilges 
as  the  ship  rolls.  These  consist  of  all  stations  for  which  the 
curvature  is  a  maximum  at  the  bilges.  They  range  from  full, 
almost  rectangular  merchant  vessel  midship  sections  to  typical 
destroyer  midship  sections. 

Note:  See  Section  2.2  for  guidelines  regarding  accuracy  of  inp\it  offsets. 
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Record  (i2a),  Horizontal  Offsets  at  Station  XA(I) 

Read  only  if  NOFF(I)  >  1  in  Record  (il). 

YA(I,J)  (NOFF(I)  reals) 

YA(I,J)  NOFF(I)  horizontal  offsets  of  station  I  (ft  or  m).  The  first  offset  point  of  each 
station  must  have  a  YA  value  of  0. 


Record  (i2b),  Vertical  Offsets  at  Station  XA(I) 

Read  only  if  NOFF(I)  >  1  in  Record  (il). 

ZA(I,J)  (NOFF(I)  reals) 

ZA(I,J)  NOFF(I)  vertical  offsets  of  station  I  (ft  or  m).  Values  are  given  as  heights  above 
the  hull  baseline. 

Notes;  1.  In  the  case  NOFF(I)  <  8  for  all  I,  the  last  point  must  be  at  the  intersection 

of  the  load  waterline  with  the  station  contour.  See  Figure  8. 

2.  In  the  case  NOFF(I)  >  8  for  some  I,  the  last  offset  point  should  be  at  the 

deck  edge.  Ideally,  the  user  should  input  8  offsets  below  the  design  waterline 
and  4  offsets  above  it.  See  Figure  9. 


Record  (i3),  Beam,  Draft  and  Area  Coefficient  at  Station  XA(I) 
Read  only  if  NOFF{I)  <  1  in  Record  (il). 

BEAM(I),  DRAFT(I),  ACOEF(I)  (3  reals) 


BEAM(I) 

DRAFT(I) 

ACOEF(I) 

Note: 


Beam  at  station  I  (ft  or  in). 

Draft  at  station  I  (ft  or  m). 

Area  coefficient  at  station  I. 

Offsets  are  generated  internally  using  either  a  Lewis  form  or.  if  the  section  is 
bulbous,  an  M.I.T.  form. 


Record  (i4),  Control  Integer  for  Load  Waterline  Calculation 
Read  only  if  NOFF(I)  >  8  for  some  I  in  Record  (il). 

IBAL  (1  integer) 

IBAL  Indicates  the  method  used  to  specify  the  loading  conditions. 

IB.4L  =  0  Input  displacement  and  LCG  (Record  (i5)). 

IB.4L  =  1  Input  midship  draft  and  trim  by  stern  (Record  (i6)). 
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Record  (i5),  Waterline  from  Displacement  and  Longitudinal  Centre  of  Gravity 
Read  only  if  NOFF(I)  >  8  for  some  I  in  Record  (il)  and  IBAL  =  0  in  Record  (i4). 
DISP,  RLCG  (2  reals) 

DISP  Ship  displacement  (tonnes  or  tons). 

RLCG  Distance  from  the  forward  perpendicular  to  the  LCG  (ft  or  m). 

Record  (i6),  Waterline  from  Midship  Draft  and  Trim  by  the  Stern 

Read  only  if  NOFF(I)  >  8  for  some  I  in  Record  (il)  and  IBAL  =  1  in  Record  (i4). 
DMID,  TRIM  (2  reals) 

DMID  Draft  at  midships  (ft  or  m). 

TRIM  Trim  by  the  stern  (positive  for  stern  down)  relative  to  bow  (ft  or  m).  This  trim  is 
the  difference  between  the  actual  trim  and  the  trim  at  the  hydrostatic  condition 
at  which  offsets  were  measured. 

Record  (j),  Positions  for  Seakeeping  Calculations  and  Slamming  Parameters 
NPOS,  THR,  SLAMEX,  (1  integer,  2  reals) 

NPOS  Number  of  positions  for  irregular  seas  seakeeping  calculations  (maximum  10). 

THR  Time  period  (hours)  for  which  number  of  slams,  slamming  pressures  and  slam¬ 
ming  forces  are  computed  for  the  stations  specified  in  Records  (jl).  If  THR  = 
0.0,  a  default  value  of  20  hours  is  used. 

SLAMEX  Exceedence  parameter  for  calculating  extreme  slamming  pressure  for  design  con¬ 
sideration.  If  SL.\MEX  =  0.0,  a  default  value  of  0.01  is  used. 

Note:  One  Record  (jl)  is  required  for  each  of  the  NPOS  positions. 

Record  (jl),  Geometry  and  Options  for  Seakeeping  Calculations 
Read  only  if  NPOS  >  0  in  Record  (j). 

XST(I),  YPOS(I),  ZPOS(I),  DECKH(l).  ISLAM(I),  IMII(I)  (NPOS  records  of  4 
reals,  2  integers) 

XST(I)  Station  number  of  position  I  where  calculations  are  to  be  done  (0  <  XST  < 
NSTOT  -  1). 

YPOS(I)  Horizontal  coordinate  (y-coordinate)  of  position  I  relative  to  the  centreline  (ft  or 
m).  If  YPOS(I)  is  not  equal  to  zero,  then  see  Note  1  for  Record  (el). 

ZPOS(I)  Vertical  coordinate  (z-coordinate)  of  position  I  relative  to  the  baseline  (ft  or  m). 

DECKH(I)  Distance  from  the  baseline  to  the  deck  edge  at  position  I  (ft  or  m).  If  the  resulting 
freeboard  is  negative,  deck  wetness  computations  are  omitted. 
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ISLAM(I)  Control  integer  for  slamming  and  deck  wetness  calculations  at  position  I. 


ISLAM  =  0 
ISLAM  =  1 

ISLAM  =  2 

ISLAM  =  3 
ISLAM  <  0 


No  slamming  or  deck  wetness  calculations  are  performed  at 
station  1.  Records  (j2)  or  (j3)  are  not  required. 

Slamming  calculations  are  performed  at  station  I  using  the 
truncated  wedge  approach.  The  form  factor  may  be  either 
given  in  Record  (j2)  or  computed  using  the  truncated  wedge 
approach  of  Reference  20.  Record  (j2)  is  required.  Deck  w’et- 
ness  calculations  are  also  performed. 

Slamming  calculations  are  performed  at  station  I  using  offsets 
given  in  Record  (j3).  Deck  wetness  calculations  are  also  per¬ 
formed. 

Deck  wetness  calculations  are  performed  at  station  I.  No  Records 
(j2)  or  (j3)  are  required. 

Propeller  emergence  calculations  are  performed  at  station  I. 
No  Records  (j2)  or  (j3)  are  required. 


IMII(I)  Control  integer  for  Mil  calculations  at  position  I. 

IMII  =  0  No  Mil  calculations  are  performed  at  position  I.  No  Record 
(jd)  is  required. 

IMII  =  1  Mil  calculations  are  performed  at  position  I.  Record  (j4)  is 
required. 

Note:  For  each  seakeeping  position,  on-  Record  (jl)  is  required  followed  by; 

ISLAM  <  0  No  Record. 

ISL.\M  =  1  Record  (j2). 

ISL.4M  =  2  Record  (j3). 

In  addition,  an  Mil  data  record  may  be  required  as  follows: 

IMII  =  0  No  Record. 

IMII  =  1  Record  (j-l). 


Record  (j2),  Geometry  for  Truncated  Wedge  Slamming  Calculations 
Read  only  if  ISLAM  =  1  in  Record  (jl). 

DEADR,  HWFB,  FFACT  (3  reals) 

DEADR  Dcadrise  angle  at  position  I  (degrees).  For  correct  slamming  force  calculations, 
DE.4DR  >  0  degrees.  If  DEADR  <  5  degrees,  the  internal  calculation  of  form 
factor  may  not  be  accurate  and  a  warning  will  be  issued.  In  this  case,  it  is 
recommended  that  ihe  form  factor  be  input  (FFACT  >  0.0)  or  that  offsets  be 
used  (ISLAM  =  2)  a.s  described  in  Record  (j3) 


HWFB  Half- width  of  flat  bottom  at  position  I  (ft  or  m). 


FFACT  Form  factor  at  position  I.  If  FFACT  <  0.0,  the  form  factor  is  computed  internally 
using  the  truncated  wedge  approach. 

Record  (j3),  Offsets  Below  l/lOth  Draft  for  Slamming  Calculations 
Read  only  if  ISLAM  =  2  in  Record  (jl). 

N,  YSTA(J)  (1  integer,  N  reals) 

N  Number  of  offsets  input  for  station  I  in  performing  slamming  calculations.  Max¬ 
imum  value  for  N  is  11. 

YSTA(J)  N-dimensional  array  whose  elements  are  the  horizontal  offsets  (distance  from  the 
centerline)  read  at  equal  vertical  intervals,  starting  at  0.1  Tx  and  continuing  down 
to  the  edge  of  the  flat  bottom  (ft  or  m).  See  Figure  10. 

Record  (j4),  Data  for  Mil  Calculations 
Read  only  if  IMIl  =  1  in  Record  (jl). 

TCLAT(i),  TCLONG(I),  TIMEOP(I)  (3  reals) 

TCLAT(I)  Lateral  tipping  coefficient  at  position  I.  If  TCLAT  <  0.  a  default  value  of  0.25  is 
used. 

TCLONG(I)  Longitudinai  tipping  coefficient  at  position  1.  If  TCLONG  <  0,  a  default  value 
of  0.17  is  used. 

TIMEOP(I)  Time  of  operation  at  position  I  (s).  If  'T'IMEOP  <  0.  a  default  value  of  60  seconds 
is  used. 

Note:  For  sliding  calculations,  the  user  should  input  TCLAT(I)  =  TCLONG(I)  =  the 
static  coefficient  of  friction. 

Record  (k),  Number  of  Bilge  Keel  Pairs 
NBKP  (1  integer) 

NBKP  Number  of  bilge  keel  pairs  (maximum  value  =  5). 

Record  (kl),  First  and  Last  Stations  Spanned  by  Bilge  Keel 
Read  only  if  NBKP  >  0  in  Record  (k). 

NFBK(I),  NLBK(I)  (NBKP  records,  2  integers) 

NFBK(I)  First  station  spanned  by  bilge  keel  I. 

NLBK(I)  Last  station  spanned  by  bilge  keel  I 
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Record  (k2),  Bilge  Keel  Offsets,  Breadths  and  Lengths 
Read  only  if  NBKP  >  0  in  Record  (k). 

YBK(J),  ZBK(J),  BKK(J),  ELBK(J)  (4  reals/record) 

YBK(J)  Horizontal  bilge  keel  offset  at  station  J  (ft  or  m). 

ZBK(J)  Vertical  bilge  keel  offset  at  station  J  (ft  or  m)  relative  to  hull  baseline. 

BKK(J)  Bilge  keel  breadth  at  station  J  (ft  or  m). 

ELBK(J)  Bilge  keel  length  at  station  J  (ft  or  m). 

Notes:  1.  One  Record  (k2)  is  required  for  each  station  spanned  by  bilge  keel  I.  That 

is,  following  each  Record  (kl),  there  must  be  (NLBK(I)  —  NFBK(I)  +  1) 
Records  (k2). 

2.  See  Figure  11  for  clarification  of  bilge  keel  inputs. 


Record  (1),  Skeg  Data 

XSK,  BSK,  ELSK  (3  reals) 

XSK  Station  number  (0  <  XSK  <  NSTOT  -  1)  of  aftmost  point  where  skeg  meets 
hull. 

BSK  Skeg  Breadth  (ft  or  m),  see  Figure  12. 

ELSK  Skeg  Length  (ft  or  m),  see  Figure  12. 

Record  (m),  Number  of  Fin  Pairs  and  Propeller  Shaft  Brackets 

NFP,  NSH  (2  integers) 

NFP  Number  of  fin  pairs  or  tank  stabilization  indicator. 

NFP  >  0  The  ship  has  NFP  fin  pairs  for  roll  stabilization  but  no  tank. 
Maximum  value  of  4. 

NFP  =  0  The  ship  has  neither  fins  nor  tank. 

NFP  <  0  The  ship  uses  a  U-tank  for  roll  stabilization  but  no  fins. 

NSH  Number  of  propeller  shaft  brackets  (maximum  value  =  4).  A  shaft  bracket  con¬ 
sists  of  2  arms,  as  show  in  Figure  13. 


Record  (ml).  Fin,  Bracket  and  Rudder  Data 

X(I),  Y(I),  Z(I),  B(I),  CR(I),  CE(I),  CLAO(I),  GAM(K),  VEFF(I)  (NFS  records  ,  9 
reals/record) 

X(I)  Station  at  whi(.  .til  I  is  located  (0  <  X  <  NSTOT  -  1). 
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Y(I) 


Z(I) 

B(I) 

CR(I) 

CE(I) 

CLAO(I) 

GAM(I) 


VEFF(I) 


Notes; 


Horizontal  coordinate  of  root  of  foil  I  (ft  or  m).  For  rudders,  if  Y(NFS)  =  0,  the 
ship  is  assumed  to  have  a  single  rudder.  If  the  rudder  stock  offset  distance  is 
positive  (Y(NFS)  >  0),  it  is  assumed  that  there  are  twin  rudders,  each  with  the 
specified  span  and  chord.  The  i,  y,  and  z  coordinates  refer  to  the  rudder  stock. 

Vertical  coordinate  of  root  of  foil  I  (ft  or  m)  relative  to  hull  baseline. 

Span  of  foil  I  (ft  or  m). 

Root  chord  of  foil  I  (ft  or  m). 

Tip  chord  of  foil  I  (ft  or  m). 

Lift  curve  slope  of  foil  I  (rad"').  If  CLAO(I)  is  unknown,  input  zero  and  the 
program  will  calculate  a  value. 

Dihedral  angle  of  foil  I  (degrees).  For  rudders,  the  dihedral  angle  is  always  taken 
as  90  degrees,  regardless  of  the  value  read  for  GAM(NFS).  Hence,  GAM(NFS)  is 
used  merely  to  specify  whether  or  not  the  rudder  is  in  the  propeii^'  E’‘nstream, 
where; 

GAM(NFS)  =  0  Rudder  not  in  propeller  slipstream,  and  the  rudder  lift  curve 
slope  is  corrected  for  boundary  layer  effects. 

GAM(NFS)  =  1  Rudder  is  in  propeller  slipstream,  and  no  boundary  layer  cor¬ 
rection  is  made. 

Correction  factor  applied  to  the  ship  forward  speed  U  in  determining  ihe  effective 
flow  velocity  V  over  foil  1.  V  =  U  x  VEFF(I).  If  foil  is  in  wake,  VEFF(I)  <  1.  If 
foil  is  in  propeller  slip-stream,  VEFF(I)  >  1.  If  VEFF(I)  <  0.0,  a  default  value 
is  taken  as  VEFF(I)  =  1.0.  For  warships,  VEFF(I)  is  typically  between  0.8  and 
1.2,  but  it  is  difficult  to  estimate  VEFF(I)  accurately.  It  is  recommended  that 
the  user  input  VEFF(I)  =  1.0  in  all  cases,  since  this  is  sufficiently  accurate  for 
engineering  purposes. 

1.  Fins,  propeller  shaft  brackets,  and  rudders  are  regarded  as  foils.  One  Record 
(ml)  is  necessary  for  each  foil  in  the  following  order; 

NFP  Records  for  the  fin  pairs  (NFP  >  0), 

NSH  records  for  the  propeller  shaft  brackets, 

1  record  for  the  rudder  or  rudders. 

That  is,  NFS  =  NFP  -|-  NSH  -i-  1  Records  (ml)  are  required. 

2.  For  NFP  fin  pairs,  coordinates  are  input  for  the  port  foil;  it  is  assumed  that 
a  counterpart  exists  on  the  starboard  side.  Note  that  y  and  z  coordinates 
apply  to  the  foil  root.  The  z  coordinate  is  input  as  height  above  baseline. 
Foil  inputs  are  illustrated  in  Figure  13. 

3.  A  shaft  bracket  consists  of  2  arms,  each  considered  as  a  foil.  For  NSH  =  1, 
input  data  only  for  the  port  arm  of  the  bracket.  For  NSH  =  2  or  4,  input 
data  for  both  arms  of  each  port  bracket,  with  the  record  for  the  outboard 
arm  preceding  that  for  the  inboard  arm.  For  NSH  =  3,  input  data  for  the 
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outboard  arm  of  tlie  port  bracket,  followed  by  the  inboard  arm  of  the  port 
bracket,  followed  by  the  port  arm  of  the  centre  bracket. 


Record  (m2),  Nominal  Fin  Angle  and  Flap  Fractional  Fin-Lift  Increase 
Read  only  if  NFP  >  0  in  Record  (m). 

BETNOM,  FLINC  (2  reals) 

BETNOM  Nominal  fin  angle  (degrees).  BETNOM  is  used  in  estimating  fin-fin  and  fin-bilge 
keel  interference  effects. 

FLINC  Fractional  increase  in  fin  lift  due  to  flap  deflection  (for  an  unflapped  fin  FLINC 

=  0.0). 

Record  (m3),  Rudder  Roll  Gains,  Natural  Frequency  and  Damping 
QFDDR,  QFDR,  QFR,  \VR,  ZETR,  WLR,  WHR  (7  reals) 

QFDDR  Rudder  roll  acceleration  gain  (sec^). 

QFDR  Rudder  roll  velocity  gain  (sec). 

QFR  Rudder  roll  gain. 

WR  Rudder  control  system  natural  frequency  (rad/sec). 

ZETR  Rudder  control  system  damping  ratio. 

WLR  Low  frequency  cut-off  for  rudder  roll  stabilizer  (rad/sec). 

WHR  High  frequency  cut-off  for  rudder  roll  stabilizer  (rad/sec). 

Note:  If  modelling  of  the  steering  system  is  not  required,  all  parameters  in  Record  (m3) 
should  be  set  to  zero. 

Record  (m4).  Rudder  Yaw  Gains 
QYDDR,  QYDR,  QYR  (3  reals) 

QYDDR  Rudder  yaw  acceleration  gain  (sec^). 

QYDR  Rudder  yaw  velocity  gain  (sec). 

QYR  Rudder  yaw  gain. 

Note:  If  modelling  of  the  steering  system  is  not  required,  all  parameters  in  Record  (m4) 
should  be  set  to  zero. 
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Record  (m5),  Fin  Roll  Gains,  Natural  Frequency  and  Damping 
Read  only  if  NFP  >  0  in  Record  (m). 

QFDD,  QFD,  QF,  \VF,  ZETF,  WHF  (6  reals) 

QFDD  Fin  roll  acceleration  gain  (sec^). 

QFD  Fin  roll  velocity  gain  (sec). 

QF  Fin  roll  gain. 

WF  Fin  control  system  natural  frequency  (rad/sec). 

ZETF  Fin  control  system  damping  ratio. 

WHF  High  frequency  cut-off  for  fin  roll  stabilizer  (rad/sec). 

Record  (m6),  Anti-Rolling  Tank  Data 
Read  only  if  NFP  <  0  in  Record  (m). 

TL,  WT,  WIT,  ZT,  ALF.4T,  HT,  RTBL.  XTST,  RHOT,  ^■S,  VSO.  WB  (12  reals) 

TL  Longitudinal  length  of  anti-rolling  tank  (ft  or  m). 

WT  Width  of  tank  connecting  duct  (ft  or  m). 

WlT  Bottom  width  of  tank  vertical  leg  (ft  or  m). 

ZT  Average  fluid  depth  in  tank  vertical  leg  (ft  or  m). 

ALFAT  Inclination  of  outside  wall  of  lank  vertical  leg  (degrees). 

HT  Height  of  tank  connecting  duct  (ft  or  m). 

RTBL  Height  above  baseline  of  bottom  of  tank  connecting  duct  (ft  or  m). 

XTST  Station  number  of  tank  location  (0  <  XTST  <  NSTOT  -  !)■ 

RHOT  Specific  gravity  of  tank  fluid. 

VS  Tank  valve  resistance  coefficient  for  frequencies  w’  >  WB. 

VSO  Tank  valve  resistance  coefficient  for  w  =  0.  Set  to  0.0  for  a  pure  passive  tank. 
WB  Break  frequency  (rad/sec).  Set  to  0.0  for  a  pure  passive  tank. 

Notes:  1.  Tank  inputs  are  illustrated  in  Figure  14. 

2.  Tank  valve  resistance  coefficient  may  be  assigned  any  value  between  zero 
and  infinity.  The  value  used  in  computations  is: 

VS  for  w  >  WB 

K  =  { 

VSO  -I-  (VS  -  VSO)..VWB  for  w-  <  W B 
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This  scheme  enables  modelling  of  a  controllable  passive  tank  in  which  damp¬ 
ing  is  effectively  increased  at  low  frequencies  in  order  to  avoid  amplification 
of  roll  by  the  stabilizer  in  quartering  seas.  For  example,  satisfactory  simula¬ 
tion  of  a  controllable  passive  tank  installed  on  a  destroyer  has  been  achieved 
by  setting  VS  =  15,  \VB  =  0.9  RNF,  and  VSO  =  25  VS. 

Record  (n),  Control  Integer  for  Next  Case 
NEXT  (1  integer) 

NEXT  Control  integer  to  end  program. 

NEXT  >  0  Input  data  for  another  case  follows  immediately  (i.e.  Records 
(a)-(n)). 

NEXT  <  0  Terminate  program  following  all  calculations  for  the  preceding 
data. 

A. 2  Changes  to  Input  Records  Since  SHIPMO  Version  4 

The  following  is  a  list  of  changes  made  to  input  records  between  SHIPMO  versions  4  and 
5.  Note  that  SHIPM05  is  capable  of  reading  older  input  files  without  modification,  as  described 
in  Section  4.3 

Hydrodynamic  End-Effect  Control  Integer  lEND  in  Record  (d) 

Control  integer  lE.N'D  is  included  in  Record  (d). 

Principal  Sea  Directions  PSDIR(I)  in  Record  (el)  for  Short-Crested  Seas 

If  responses  in  short-crested  seas  are  requested  (NSD  =  0  or  .4NGLE  >  0  in  Record  (e)). 
attention  must  be  given  to  PSDIR(I).  If  all  seaheeping  positions  arc  on  the  centreline  (VPOS(I) 
=  0  for  all  I  in  Record  (jl)).  PSDIR  must  cover  the  range  0-180  degrees.  If  any  seakeeping 
position  is  off  centreline,  PSDIR  must  cover  the  range  from  0  to  a  minimum  of  315  degrees. 
PSDIR  can  be  in  ascending  or  descending  order,  with  a  maximum  increment  of  45  degrees 
between  values.  SHIPM05  checks  to  see  if  these  criteria  are  satisfied,  and  terminates  execution 
if  the  input  sea  directions  need  to  be  modified. 

Centre  of  Gravity  Location  KG  in  Record  (h) 

The  height  of  the  CG  is  given  relative  to  the  ba.seline. 

Y  Offset  of  Seakeeping  Position  YPOS(I)  in  Record  (jl) 

The  y  offset  YPOS(I)  for  a  scakeeping  position  is  included  in  Record  (jl ). 


Vertical  Coordinate  of  Seakeeping  Position  ZPOS(I)  in  Record  (jl) 

The  vertical  coordinate  ZPOS(I)  is  be  given  relative  to  the  ship  baseline. 

Vertical  Coordinate  of  Deck  Edge  DECKH(I)  in  Record  (jl) 

The  vertical  coordinate  of  the  deck  edge  DECKH(I)  is  given  relative  to  the  ship  baseline. 

Motion-Induced  Interruption  Control  Integer  MII(I)  in  Record  (jl) 

Control  integer  MII(I)  is  included  in  Record  (jl). 

Height  of  Bottom  of  Anti-Roll  Tank  RTBL  in  Record  (m6) 

The  height  of  the  bottom  of  the  tank  connecting  duct  RTBL  is  given  relative  to  the  ship 
baseline. 

Longitudinal  Location  of  Anti- Roll  Tank  XTST  in  Record  (m6) 

The  longitudinal  location  of  the  anti-roll  tank  is  input  as  station  number  XTSl . 
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A. 3  Sample  Input  File 


Record 

(e) 

FRIGATE 

EXAMPLE 

< - Title 

Record 

(b) 

0  10  1 

2 

2  0 

<— 

- Control  iii-egers 

Record 

(c) 

.2  2.0 

0.1 

< - 

Vav*  frequencies 

Record 

(d) 

0  0  1  5 

.8  0. 

2  0 

0  < 

Encounter  frequencies.  Method,  l£ID 

Record 

(e) 

7.  90. C 

i  < 

Roaber 

of  sea  directions,  ASCLE 

Record 

(el) 

0  30  60  90  120  ISO  180.0 

<— 

See  directions 

Record 

(f) 

1  <— 

- loaber 

of  sesssjrs 

Record 

(f2) 

10.66, 

10.0 

<-  — 

—  Sig&ificsnt  oaee  height,  period 

Record 

(8) 

18  18 

<- 

-  Ship  speeds 

Record 

(h) 

356.0  18.09 

2.09  C 

>.0  0.0 

0.35 

0.25  0  < - Basic  ship  data 

Record 

(i) 

21  20 

.  1. 

< 

loaber  of  stations,  scaling  factors 

Record 

(il) 

1  12 

1 

-  StstioD  Bimber  1  data 

Record 

(i2e} 

0.0 

0.46 

1.19 

1 .82 

2.33 

2.57  2.84  3.50  4.50  6.55  9.25  12.79 

Record 

(i2b) 

0.0 

0.00 

2.77 

6.60 

9.93 

11.21  12.61  15.53  18.84  23.84  28.78  34.01 

Record 

(il) 

2  12 

1 

Record 

(i2e) 

0.0 

0.49 

1.60 

2.76 

3.61 

4.49  5.04  6.56  8.32  11.19  13.65  15.57 

Record 

(i2b) 

0.0 

0.0 

1.44 

3.62 

5.96 

8.70  10.56  15.05  19.40  25.42  29.81  32.65 

Record 

(il) 

3  12 

1 

Record 

(i2a) 

0.0 

0.49 

3.07 

4.55 

5.77 

6.90  8.24  9.58  11.19  13.40  15.19  17.16 

Record 

(i2b) 

0.0 

0.0 

2.01 

3.69 

6.06 

8.49  11.76  14.97  18.87  23.68  27.30  31.40 

Record 

(il) 

4  12 

1 

Record 

(i2a) 

0.0 

0.49 

3.79 

6.43 

8.75 

10.19  11.53  12  72  14  00  15  38  16.92  18  09 

Record 

(i2b) 

0.0 

0.0 

1.54 

3.86 

7.13 

9.97  12.77  15.58  18.86  22  32  26.88  30.34 

Record 

(il) 

5  12 

1 

Record 

(i2a) 

0.0 

0.55 

4.59 

8.68 

11.01 

12.75  14.52  16.01  16.37  17.33  17.92  18  93 

Record 

(i2b) 

0.0 

0.0 

1 .30 

4.18 

6.86 

9.76  13.47  17.34  18  82  22.12  25.01  29.30 

Record 

(il) 

6  12 

1 

Record 

(i2a) 

0.0 

0.52 

7.18 

11.60 

14.70 

15.80  16.62  17.90  18.25  18  81  19.24  19.65 

Record 

(i2b) 

0.0 

0.0 

1 .93 

4.94 

8.78  10.59  12.55  16  66  18.65  21.67  25.21  28.14 

Record 

(il) 

7  12 

0 

Record 

(i2a) 

0.0 

0.52 

8.78 

12.59 

15.29 

17.10  18.14  19.11  19.64  20.06  20.22  20,27 

Record 

(i2b) 

0.0 

0.0 

2.10 

4.01 

6.39 

9.04  11.22  14.52  18.60  23.83  26.69  27.53 

Record 

(il) 

8  12 

0 

Record 

(i2a) 

0.0 

0.49 

8.87 

13.40 

16.23 

18.61  19.71  20.40  20.47  20.65  20.66  20.74 

Record 

(i2b) 

0.0 

0.0 

1  .91 

3.46 

5.27 

8.46  11.66  15.12  18.49  22  64  25  47  27  16 

Record 

(il) 

9  12 

3 

Record 

(i2a) 

0.0 

0.49 

8.90 

14.66 

17.43 

19.36  20.39  20.84  20.99  20.94  20.92  20.94 

Record 

(i2b) 

0.0 

0.0 

1 .91 

3.49 

5.16 

7.53  10.58  13.75  19.04  22.18  24.92  27.04 

Record 

(il) 

10  12 

3 

Record 

(i2a) 

0.0 

0.46 

12.67 

16.51 

18.66 

19.73  20.51  20.87  21.00  20.97  20.92  20.94 

Record 

(i2b) 

0  0 

0.0 

2.65 

4.16 

5.98 

7.69  10.03  12.98  16.12  20.72  24.06  26.50 

Record 

(il) 

11  12 

3 

Record 

(i2a) 

0  0 

0.43 

12.67 

16  39 

18.53 

19.54  20.36  20  77  20.97  20  94  20,92  20.94 

Record 

(i2b) 

0  0 

0.0 

2.71 

4.32 

6.16 

7.66  10.15  13.01  16.12  20.72  24.06  26.50 

Record 

(il) 

12  12 

3 

Record 

(i2a) 

0.0 

0.43 

9  90 

14.33 

17.40 

18.87  20.08  20.62  20  91  20.91  20.92  20  91 

Record 

(i2b) 

0.0 

0.0 

2  14 

3.65 

5.66 

7.28  9.90  12.95  16.12  20.72  24  06  26.5 

Record 

(il) 

13  12 

3 

Record 

(i2a) 

0.0 

0.43 

11.53 

14  72 

16  69 

16.74  19.68  20.41  20.81  20.85  20.86  20,85 

Record 

(i2b) 

0.0 

0.0 

2  99 

4.34 

5.66 

7.90  9.87  12.64  16.12  20.72  24  09  26.50 

K«cord 

(il) 

14  12 

3 

Record 

(i2a) 

0.0 

0.40 

9  56 

13.95 

17.48 

18.76  19.68  20.28  20.60  20.76  20.77  20  75 

E«cord 

(i2b) 

0  22 

0.22  2.82  4.74  7.34  9.05  11.14  13.63  16.12  20.72  24.03  26.50 

Record 

(il) 

15  12 

2 

Record 

(i2a) 

0.0 

0.40 

8  36 

14  56 

16  95 

18.51  19.39  20.02  20.26  20  48  20.46  20.45 

E«cord 

(i2b) 

0.96 

0.96  3.40  6.17 

7.99  10  14  12.41  14.93  16  15  20.72  24  06  26  5 

Record 

(il) 

16  12 

2 

Record 

(i2a) 

0.0 

0.51 

6  41 

11.21 

14.58 

16.45  17.89  19  16  19.40  19.68  19  75  19  80 

Record 

(i2b) 

2.70 

2.70  4.93  5.97  7.66  9.26  11.35  15.27  16.08  20.71  24.05  26.5 

Record 

(il) 

17  12 

2 

Record 

(i2a) 

0.00 

0.46  7.66  12.55  14.85  16.35  17.02  17.41  18  10  18  60  18.77  18  84 

Record 

(i2b) 

5.41 

5.41 

7.06  8.55 

9.99  11.56  12  86  13.86  16  08  20  64  24.02  26  5 

Record  (il) 
Rocord  (i2«) 
Record  (13b) 
Record  (11) 
Record  (13e) 
Record  (13b) 
Record  (11) 
Record  (13a) 
Record  (13b) 
Record  (14) 
Record  (16) 
Record  (j) 
Record  (jl) 
Record  (j3) 
Record  (j4) 
Record  (k) 
Record  (kl) 
Record  (k3) 
Record  (k3) 
Record  (k3) 
Record  (k2) 
Record  (k3) 
Record  (1) 
Record  (■) 
Record  (■!) 
Record  (ml) 
Record  (■!) 
Record  (k3) 
Record  (>4) 
Record  (n) 


18  13  2 

0.0  0.38  10.30  12.36  13. SS  14.62  15.38  15.78  16.33  17.00  17.38  17.64 

8.48  8.48  10.06  10.73  11.48  13.42  13.63  14.69  16.10  20.66  24.01  26.5 

19  12  0 

0.0  6.56  8.93  10.23  11.33  11.97  12.61  13.25  13.82  14.88  15.59  16.10 

10.66  10.75  10.89  11.37  12.02  12.65  13.59  14.83  16.08  20.65  24.00  26.5 
30  12  0 

0.03  4.74  6  25  7.23  8.15  8.80  9.65  10.17  10.78  12.20  13.29  14.13 
10.66  10.81  10.97  11.32  11.85  13.45  13.63  14.75  16.06  30.63  24.02  36.5 

1  < -  Control  integer  for  load  naterline 

13.69  0.5  < - Draft  at  aidakips,  trin  bj  stern 

1  20.0  .01  < - -  Inaber  of  positions 

3.0  0.0  29.27  33.29  1  1  < -  Position  for  seakeeping  calcnlation 

34.5  .5  0.  < -  Data  for  slasaing  calcnlation 

0.25  0.17  60.0  < -  Tipping  coefficient,  tiae  of  operation 

1  < -  Inaber  of  bilge  keel  pairs 

10  14  < -  First  and  last  stations  spanned  bj  bilge  keel 

18.960  7.010  3.000  12.100  < -  Bilge  keel  offsets 

18.780  6.100  2.000  17.800 
18.470  6.000  3.000  17.800 
18.350  6.380  2.000  17.800 
18.210  7.130  3.000  21.900 
17.87  3.900  60.000  < -  Skeg  data 

0  3  < -  Inaber  of  fin  pairs,  propeller  shaft  brackets 

19.30  10.200  10.500  7.400  2.625  2.625  0.000  79.000  1.000  <—  Outboard  ara 
19.30  3.800  9.700  7.350  2.625  2.625  0.000  65.000  1.000  <—  Inboard  ara 

19.16  6.500  5.850  9.920  7.5  6.017  0.0  1.0  1.000  <—  Rudder  data 

0.000  0.000  0.000  0.000  0.000  0.000  0.000  < -  Rudder  roll  gains 

0.000  0.000  0.000  < Rudder  Jae  gains 

0 


41 


A. 4  Input  Format  for  Hindcast  Directional  Spectrum 

When  a  hindcast  directional  spectrum  is  selected  by  the  user  (ISPEC  =  G  in  Record  (b)). 
spectral  data  is  read  from  a  formatted  file  specified  in  Record  (fG).  The  contents  of  the  file  are 
as  follows: 

Record  1  -  Labels  for  File  Descriptors 
Format  (A89) 

LABELS  Characterise  Labels  for  values  on  following  line. 


Record  2  -  File  Descriptors 

Format  (12,  12,  IX,  12,  12,  IX,  16.  IX,  F8.3,  F8.2,  2X,  F6.2.  2X.  F6.2.  2X.  F6.3, 
4F8.2) 


lYR 

Integer 

Last  two  digits  of  year. 

IMON 

Integer 

Month. 

IDAY 

Integer 

Day  of  month. 

IHR 

Integer 

Hour  (GMT). 

IGRP 

Integer 

Grid  point  number. 

LAT 

Real 

Latitude  (degrees). 

LONG 

Real 

Longitude  (degrees). 

WSPD 

Real 

Wind  speed  (m/s). 

WDIR 

Real 

Wind  direction  (degrees)  from  which  wind  is  approaching 

USHEAR 

Real 

Wind  shear  velocity  (m/s). 

HS 

Real 

Significant  wave  height  (m). 

TS 

Real 

Significant  wave  period  (m). 

VMD 

Real 

Mean  direction  (to)  of  energ.v  propagation  (degrees). 

TP 

Real 

Peak  wave  period  (s). 

Record  3  - 

Blank 
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Record  4  -  Wave  Frequencies 

Format  (A9,  3X,  15F8.5) 

LABEL  Character^g  Line  label  “FREQUENCY'^ 

FREQ(15)  Real  Array  of  circular  frequencies  for  energy  components  (Hz). 

Record  5  -  Direction  Label 

Format  (A9) 

LABEL  Character*9  Line  label  “DIRECTION". 

Remaining  Records  -  Directional  Spectral  Energies  (Maximum  24  Records) 

Format  (F7.1,  “iX,  15F8..5) 

DIR  Real  Compass  direction  (to)  for  energy  components  on  line  (degrees). 

E(15)  Real  Total  energy  for  given  direction  and  frequency  (m^). 

Note:  Each  spectral  energy  record  covers  a  directional  range  of  15  degrees:  however,  no 

records  are  given  for  directions  with  very  low  energies. 


A. 5  Sample  Input  Hindcast  Directional  Spectrum 
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A. 6  Input  Format  for  User  Input  Directional  Spectrum 


When  a  user  input  directional  spectrum  is  selected  (ISPEC  =  7  in  Record  (b)).  spectral 
data  is  read  from  a  formatted  file  specified  in  Record  (f7).  The  contents  of  the  file  are  as  follows: 

Record  1  -  Spectrum  Title 

SPTITLE  (Columns  1-80) 

SPTITLE  Alphanumeric  title  (maximum  of  80  characters)  which  is  written  on  output. 

Record  2  -  Spectrum  Parameters 

ISEAIN,  ITOP,  HS,  TP  (2  integers,  2  reals) 

ISEAIN  System  of  units  for  significant  wave  height  and  spectral  density: 

ISEAIN  =  0  British  units 
ISEAIN  =  1  Metric  units 

The  units  of  wave  height  (ft  or  m)  and  spectral  density  (ft^/(rad/s)/rad  or 
m^/(rad/s)/rad)  may  be  chosen  independently  of  the  value  IN’  specified  in  Record 
(b)  of  the  main  input  file. 

ITOP  Number  of  frequencies  for  spectral  density  input  (maximum  30). 

HS  Significant  wave  heiglit  (ft  or  m). 

TP  Characteristic  wave  period  (ft  or  m)- 


Record  3  -  Spectrum  Frequencies 
AW(I)  (ITOP  reals) 

AW(I)  Frequencies  for  spectral  energy  densities  (rad/s). 

Records  4-39  -  Sea  Direction  and  Spectral  Energy  Densities 

DIR,  (ENSP(IDIR,  J),  J  =  I,  ITOP)  (36  records,  1-l-ITOP  reals/record) 

DIR  Sea  direction  (from)  for  spectral  energy  densities  (degrees). 

ENSP(IDIR.  J)  Spectral  energy  densities  (m^/(rad/s)/rad  or  ft^/(rad/s)/rad). 

Note:  Sea  directions  DIR  must  begin  at  0  degrees  and  increase  in  10  degree  increments 
to  a  maximum  of  350  degrees.  The  program  will  terminate  if  this  convention  is 
not  followed. 
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B  Program  Output 

This  Appendix  gives  a  detailed  description  of  line  printer  and  post-processor  output  from 
SHIPMO. 


B.l  Printer  Output  Description 

Appendix  B.3  contains  the  program  output  which  results  from  the  sample  input  of  Ap¬ 
pendix  A. 3.  Since  the  output  is  self-explanatory,  only  a  brief  description  is  required. 

The  first  eight  pages  are  devoted  to  listing  input  data  and  results  of  hydrostatic  calculations. 
This  is  necessary  for  checking  purposes.  On  the  first  page,  information  is  given  about  Records 
(a)  to  (h).  If  a  measured  or  hindcast  spectrum  was  used,  the  spectral  data  would  have  been 
listed  on  page  two  of  the  output.  Pages  two  and  three  list  hull  offsets  and  the  fourth  and  fifth 
pages  are  line  printer  plots  of  the  body  plan.  These  provide  a  valuable  visual  check  of  the 
offset  data.  The  program  always  outputs  the  offsets  up  to  the  load  waterline,  even  when  the 
input  file  contains  offsets  up  to  the  deck  edge.  The  output  offsets  are  also  adjusted  to  account 
for  any  trim.  Note  that  the  offset  data  arc  not  output  when  ISAVE  =  2.  The  sixth  page 
contains  hydrostatic  quantities  computed  from  the  offsets,  as  well  as  data  on  specified  positions 
for  seakeeping  calculations.  On  the  seventh  page,  bilge  keel.  skeg.  fm.  .\-bracket.  rudder,  and 
control  inputs  are  listed.  .4nti-rolling  tank  inputs  would  also  be  listed  on  this  page,  if  applicable, 
as  well  as  certain  tank  parameters  such  as  tank  natural  frequency.  .41so  on  this  page  arc  the  roll 
radius  of  gyration  and  roll  natural  frequency  (rad/s). 

The  next  fourteen  pages  contain  the  motions  and  roll  damping  coefficients  in  regular  waves 
at  the  seven  sea  directions  specified.  To  minimize  unnecessary  output,  regular  wave  responses 
for  headings  greater  than  180  degrees  arc  not  output  because  of  symmetry.  For  each  heading, 
the  motion  responses  are  followed  by  roll  damping  coefficients  if  they  have  been  requested.  On 
the  output,  \V  is  wave  frequency  and  WE  is  frequency  of  encounter,  both  in  rad/s.  W.L./L  is 
wavelength  divided  by  ship  length.  Phases  are  expressed  relative  to  wave  elevation  at  the  CG: 
e.g.,  sway  phase  =  90  degrees  means  sway  achieves  it.c  maximum  positive  value  at  a  point  on 
the  wave  cycle  90  degrees  before  wave  elevation  at  the  CG  achieves  its  maximum  positive  value. 
Note  that  SHIPMO  phase  angles  can  be  converted  to  the  ITTC  convention  [7]  as  follows: 

ITTC  SHIPMO 

Surge  phase  =  Surge  phase  -I-  180  degrees 

Sway  phase  =  Sway  phase 

Heave  pha.se  =  Heave  pha.se 

Roll  pha.se  =  Roll  phase  180  degrees 

Pitch  phase  =  Pitch  phase 

4'a\\  phase  -  4’aw  phase 
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ROLANG  is  the  roll  angle  at  whicli  tlip  roll  damping  coefficients  were  computed.  If  tlie 
program  is  computing  the  motions  in  irregular  waves,  then  for  each  sea  state  and  heading. 
ROL.ANG  is  the  average  roll  amplitude,  wliich  equals  1.25  times  the  RMS  roll.  If  only  regular 
wave  computations  are  being  performed,  then  for  each  heading  and  wave  frequency.  ROL.ANG 
equals  the  wave  slope  (specified  in  Record  (fl))  times  the  roll  response  per  unit  wave  slope. 

The  translational  and  angular  amplitudes  of  the  motions  are  non-dimensionaJized  by  wave 
amplitude  and  maximum  wave  slope,  respectively. 

The  output  of  motions  and  roll  damping  coefficients  in  regular  w’aves  is  optional,  and  may  be 
deleted  by  setting  IRESP  =  0  in  Record  (b).  This  option  would  be  exercised  when  only  motions 
in  irregular  seas  are  of  interest,  as  is  often  the  case.  Note  further  that  although  frequency 
response  is  computed  for  each  seaway  specified  in  Records  (f2),  vf4),  or  (f5),  frequency  response 
is  output  only  for  the  first  seaway  when  NSE.A  >  1.  This  eliminates  unnecessary  output,  since 
frequency  responses  computed  for  different  seaways  would  have  noticeable  differences  only  in  roll, 
where  peak  response  is  inversely  proportional  to  roll  damping  which  is  directly  proportional  to 
root  mean  square  (RMS)  roll.  The  reader  is  reminded  that  in  the  case  of  a  measured  seaway 
spectrum  (ISPEC  =  2),  hindcast  spectrum  (ISPEC  =  6).  or  user  input  directional  spectrum 
(ISPEC  =  7),  NSE.A  must  equal  1. 

The  remainder  of  the  output  presents  the  results  of  irregular  sea  calculations.  In  addition  to 
the  R.MS  responses,  zero-crossing  periods  are  now  output.  If  a  non-directional  spectrum  has  been 
selected  (quadratic  regression.  Bretschneider.  JONS\V.-\P.  measured  point,  or  Ochi  and  Hubble  6 
parameter,  ISPEC  =  0-1).  seakoeping  data  are  given  for  unidirectional  and  short-crested  seas  at 
each  of  the  principal  sea  directions  (lieadings)  specified  in  Rccoid  (el).  If  symmetry  exists  (i.e. 
a  point  lies  on  the  ship  centreline),  output  is  not  given  for  headings  greater  than  ISO  degrees. 
For  the  directional  spectra  (10  parameter,  hindcast.  or  user  input  directional.  ISPEC  =  5,  6  or 
7).  the  motion  responses  are  not  symmetric  even  for  positions  on  the  centreline,  and  seakeeping 
data  are  given  only  for  short-crested  seas  for  default  ship  compass  headings  of  0-3C0  degrees. 

The  output  RMS  values  begin  with  the  ship  translational  and  rotational  motions  at  the 
CG,  as  well  as  RMS  rudder,  fin,  or  tank  angles  due  to  stabilizer  activity.  Next,  the  results  of 
seakeeping  calculations  for  the  positions  specified  in  Record  (jl)  are  listed.  \’ertical,  lateral, 
and  longitudinal  motions  are  given  for  all  seakeeping  positions.  Motion-induced  interruption 
(or  sliding)  statistics  and  forces  relative  to  loc;:l  ship  axes  are  gi\en  for  positions  with  MII(I) 
=  1.  Finally  deck  wetness  and  slamming  computations  (ISk.^M)!)  >  1 )  or  propeller  emergence 
computations  (ISL.AM(I)  <  0).  and  or  deck  wetness  calculations  (freeboard  >  0)  may  be  given 
for  seakeeping  positions.  Slamming  computations  are  performed  as  described  in  Reference  20. 
The  output  is  as  follows: 

•  RMS  relative  vertical  displacement  and  velocity. 

•  vibration  ride  quality  index  (see  Reference  19). 

•  probability  of  deck  wetness  at  any  instant. 

•  expected  number  of  deck  wetnesses  per  hour. 

•  probability  of  keel  emergence  at  any  instant. 


•  expected  keel  emergences  per  hour, 

•  most  probable  maximum  slam  pressure  during  TIIR  hours. 

•  slam  pressure  with  SLAMEX  probability  of  exceedence  during  THR  hours, 

•  most  probable  maximum  slam  force/unit  length  during  THR  hours, 

•  slam  force/unit  length  with  SLAMEX  probability  of  exceedence  during  THR  hours. 

The  following  output  values  are  given  for  propeller  emergence  computations: 

•  RMS  relative  vertical  displacement  and  velocity, 

•  probability  of  propeller  emergence  at  any  instant, 

•  expected  number  of  propeller  emergences  per  hour. 

In  the  regular  sea  case,  the  responses  are  output  in  nondimensional  form,  but  in  the  irregular 
sea  case  most  of  the  output  is  dimensional.  Hence,  when  running  SHIPM0.5  with  model  data, 
care  must  be  taken  to  interpret  the  irregular  sea  results  correctly. 

B.2  Post-Processor  Output 

In  program  SHIPMO.  frequency  responses  and  RMS  motions  may  be  stored  on  disk  for 
post-processing.  This  is  accomplished  by  specifying  IPPF  =  1  in  Record  (b).  and  giving  the 
disk  file  a  name  in  Record  (bl).  Data  will  then  be  written  into  a  sequential,  unformatted  binary 
file. 

Contents  of  the  post-processing  file  arc  given  below.  Variables  with  numbers  in  parentheses 
beside  them  (e.g.  TITLE(IO))  represent  arrays. 

Record  1  -  General  Specifications  from  Program  Input 


TITLE(IO) 

Char*10 

.Alphanumeric  title  with  time  and  date  appended 

NSP 

Integer 

Number  of  ship  speeds. 

UKMIX 

Real 

Lowest  ship  speed  (knots). 

DUK 

Retd 

Ship  speed  increment  (knots). 

NSEA 

Integer 

Number  of  seaways. 

ISEA(IO) 

Integer 

N.ATO  sea-states  for  seaways. 

HSW(IO) 

Real 

Significant  wave  heights  (ft  or  in). 

TS\V(10) 

Real 

Representative  periods  (s). 

NSD 

Integer 

Number  of  sea  directions 

-10 


PSDIR(36) 

Rea] 

Sea  directions. 

ANGLE 

Rea] 

Spreading  angle  (degrees). 

ISPEC 

Integer 

Spectrum  control  integer. 

Record  2  - 

General  Specifications  (Continued) 

NW 

Integer 

Number  of  wave  frequencies. 

WMIN 

Rea] 

Lowest  wave  frequency  (rad/s). 

DW 

Rea] 

Wave  frequency  increment  (rad/s). 

10  UT 

Integer 

Output  units  control  integer  (lOUT  =  0  for  FPS,  lOUT  =  1  for 
SI).  lOUT  determines  units  of  all  post-processing  output. 

Record  3  - 

Data  for  Seakeeping  Positions 

NPOS 

Integer 

Number  of  positioins. 

XST(IO) 

Rea] 

Position  station  numbers. 

YPOS(IO) 

Real 

Position  y-coordinates  (ft  or  m). 

ZPOS(IO) 

Rea] 

Position  z-coordinates  relative  to  CG  (ft  or  m). 

DRST(IO) 

Real 

Draft  at  positions  (ft  or  m). 

FORM(IO) 

Real 

Position  slamming  form  factors.  The  slamming  pressure  at  a  po¬ 
sition  is  obtained  as  follows: 

P(I)  =  0.5px  FORM(I)  X 

where  Ur  is  the  relative  velocity  (ft/s  or  m/s)  between  tlie  water 
surface  and  the  keel. 

FREEB(IO) 

Real 

Freeljoard  witliout  correction  for  dynamic  waterline  (ft  or  m). 

SLFF(IO) 

Rea] 

Dimensional  position  slamming  force  factor  (ft  or  m).  The  sec¬ 
tional  force  (i.e.  force/unit  length)  is  obtained  from  the  slamming 
pressure  as  follows: 

F(I)  =  SLFF{I)  X  P{J) 

SLAMEX 

Real 

Exceedence  probability  for  extreme  slamming  pressures  and  forces. 

THR 

Rea] 

Time  period  for  slamming  (hours). 

The  remaining  refolds  are  reiicated  for  each  of  tlie  NSP  ship  speeds. 
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Next  NSD  Records  -  Regular  Response  Data  for  NSD  Sea  Directions 

AMPHAS(40,  8)  Complex  Regular  wave  response  of  CG.  The  real  component  is  the 

dimensionless  amplitude  (normalized  by  wave  amplitude  or 
maximum  wave  slope)  and  the  imaginary  component  is  the 
phase  angle  (degrees).  Array  components  are  AMPHAS(I\V, 
JMODE),  where  IW  is  the  frequency  and  JMODE  is  the 
mode.  The  8  modes  are  as  follows: 


1.  Surge 

2.  Sway 

3.  Heave 

4.  Roll 


5.  Pitch 

6.  Yaw 

7.  Rudder 

8.  Stabilizer 


Note:  If  there  is  more  than  1  seaway,  AMPHAS  is  written  only  for  the  first  seaway. 


Next  NSEA  Records  -  RMS  Motions  for  NSEA  Seaways 

RMS(36,8,3)  Real  RMS  motions  at  the  CG.  .4rray  components  are  R.MS(JSD, 

KMODE.  LD\’A),  JSD  is  the  sea  direction.  KMODE  is  the 
mode  as  given  for  .4MP1I.\S.  and  LD\’.A  is  the  type  of  mo¬ 
tion.  The  3  motion  types  are  as  follows: 

1.  Displacement  (ft  or  m) 

2.  \'elocity  (ft/s  or  m/s) 

3.  Acceleration  (g's) 

RMSX(36,10,22)  Real  RMS  motions  at  seakeeping  positions.  .4rray  components 

are  RMS(JSD.  KPOS,  LMODE).  where  KPOS  is  the  sea¬ 
keeping  position  and  LMODE  is  the  mode.  The  modes  for 
seakeeping  positions  are  as  follows: 

1.  4’ertical  displacement  (ft  or  m) 

2.  \'ertical  velocity  (ft/s  or  m/s) 

3.  Vortical  acceleration  (g’s) 

4.  Relative  displacement  (ft  or  m) 

5.  Relative  velocity  (ft/s  or  m/s) 

G.  \’ibration  ride  quality  index 

7.  Lateral  displacement  (ft  or  m) 

8.  Lateral  velocity  (ft/s  or  m/s) 

9.  Lateral  acceleration  (g's) 

10.  Lateral  force  estimator  (g's) 

11.  Port  tipping  estimator  function  (g’s) 

12.  Time  derivative  of  11  (g's/s) 

13.  Starboard  tipjring  estimator  function  (g's) 
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14.  Time  derivative  of  13  (g's/s) 

15.  Longitudinal  displacement  (ft  or  m) 

16.  Longitudinal  velocity  (ft/s  or  m/s) 

17.  Longitudinal  acceleration  (g’s) 

18.  Longitudinal  force  estimator  (g’s) 

19.  Aft  tipping  estimator  function  (g’s) 

20.  Time  derivative  of  19  (g’s/s) 

21.  Forward  tipping  estimator  function  (g’s) 

22.  Time  derivative  of  21  (g’s/s) 

Note:  For  unidirectional  spectra  (ISPEC  =  0-4),  RMS  and  RMSX  are  only  writ¬ 
ten  for  long-crested  seas.  For  the  10  parameter  and  hindcast  spectra  (IS¬ 
PEC  =  5,6),  RMS  and  RMSX  are  written  only  for  directional  seas,  with 
PSDIR(JSD)  being  the  ship  compass  heading  for  RMS  and  RMSX. 

In  summary,  the  post-processing  file  contains  a  total  of  3  -f  NSPx(XSD  -)-  NSEA)  records. 


B.3  Sample  Printer  Output 


SBIPROS  —  SHIP  KOTIOIS  II  OBLIQUE  SEiS  PROGRiH  VERSIOI  S  O,  DECEMBER  1991 

FRIGRTE  EUHPLE  < - Titl#  10:43:14  3-iPR-92 


DIHEISIOIS  REQUESTED  :  IIPUT  -  FPS  OUTPUT  -  SI 
ROTIOIS  COMPUTED  FOR  :  <  SiLT>  HITEH 
SEiVlY  SPECTRUM  :  BRETSCBIEIDER 

DTIRMIC  SHELL-UP  CORRECTIOI  :  HO 
HAVE  PROFILE  CORRECTIOI  :  YES 

REGULAR  RESPOISE  PRIIT-OUT  :  YES 
ROLL  DAMPIIG  PRIIT-OUT  :  YES 
OUTPUT  DATA  SAVED  01  FILE  :  10 

HAVE  FREQUEICY  RAIGE  FOR  REGULAR  RESPOISE  :  0.300  TO  3.000  II  STEPS  OF  0.100 

BULL  SECTIOIAL  ADDED  MASS  A  DAMPIIG  CONPUTATIOIS 
EICOUITER  FREQUEICY  RAIGE  :  0.100  TO  S.900  II  STEPS  OF  0.300 

METHOD  USED  :  CLOSE-FIT 
EID  EFFECT  BYDRODYIAMIC  TERMS  USED  :  10 

7  SEA  DIRECTIOIS  FOR  MOTIOI  CALCULATIQIS 

0.0  30.0  60.0  90.0  130.0  ISO.O  180.0 

SPREADIIG  AIGLE  FOR  SHORT -CRESTED  SEA  SPECTRUM  «  90.00 

■0  OF  SCAHAYS  COISIOEREO  >  1 

SIG  HAVE  IT  HAVE  PERIOD  SEA  STATE 
FT  SEC 

10.6600  10.0000  5 

1  SHIP  SPEEDS  :  18.00  TO  18.00  II  STEPS  OF  1.00  HOTS 


SHIP  LEIGTH  «  356.00  FT 
IT  OF  CG  ABOVE  BASELIIE  >  18.09  FT 

METACEITRIC  IT  »  2 . 09  FT 

ROLL  lATURAL  FREQUEICY  >0.000 
ROLL  RADIUS  OF  GYRATIOI  /  BEAM  >  0.350 
YAH  RADIUS  OF  GYRATIOI  /  LHGTB  «  0.350 

IRG  >  0 

MIDSHIPS  DRAFT  >  13.690  FT 
TRIM  BY  THE  STERI  »  0.500  FT 
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SHIP  OFFSETS  (I*  FT)  AT  21  STiTIOfS 
SCILIIG  FACTORS  :  Y  «>  1.000  Z  »>  1.000 

TRIMKED  OFFSETS  FROM  HIDSHIPS  DRAFT  AID  TRIM 


STATIOI 

0 

I 

=  0.00 

lEDDY 

«  0 

Y- 

0.000 

Z- 

13.690 

STATIOI 

1 

X 

»  17.80 

lEDDY 

»  1 

Y- 

0.000 

0.460 

1.190 

1.820 

2.330 

2.670 

2.840 

3.033 

Z- 

0.225 

0.225 

2.995 

6.825 

10.155 

11.435 

12.835 

13.690 

STATIOl 

2 

X 

•  35.60 

lEDDY 

«  1 

Y- 

0.000 

0.490 

1.600 

2.760 

3.610 

4.490 

5.040 

6.032 

Z- 

0.200 

0.200 

1.640 

3.820 

6.160 

8.900 

10.760 

13.690 

STATIOI 

3 

X 

»  53.40 

lEDDY 

«  1 

Y- 

0.000 

0.490 

3.070 

4.550 

5.770 

6.900 

8.240 

8.973 

Z- 

0.17S 

0.175 

2.185 

3.865 

6.235 

8.665 

11.935 

13.690 

STATIOI 

4 

X 

-  71.20 

TEDDY 

«  1 

Y- 

0.000 

0.490 

3.790 

6.430 

8.750 

10.190 

11.530 

11 .856 

Z- 

O.ISO 

0.150 

1 .690 

4.010 

7.280 

10.120 

12.920 

13.690 

STATIQI 

5 

I 

-  89.00 

lEDDY 

*  1 

Y- 

0.000 

0.550 

4.590 

8.680 

11.010 

12.750 

14.520 

14.557 

Z- 

0.12S 

0.125 

1.425 

4.305 

6.985 

9.885 

13.595 

13.690 

STATIOI 

6 

X 

•  106.80 

lEDDY 

«  1 

Y- 

0.000 

0.520 

7.180 

11.600 

14.700 

15.800 

16.620 

16.944 

Z- 

0.100 

0.100 

2.030 

5.040 

8.880 

10.690 

12.650 

13.690 

STATIQI 

7 

X 

«  124.60 

lEDDY 

■  0 

Y- 

0.000 

0.520 

8.780 

12.590 

15.290 

17.100 

18.140 

18.844 

Z- 

0.075 

0.075 

2.175 

4.085 

6.465 

9.115 

11 .295 

13.690 

STATIOI 

8 

X 

»  142.40 

lEDDY 

«  0 

Y- 

0.000 

0,490 

8.870 

13.400 

16.230 

18.610 

19.710 

20.105 

Z- 

0.050 

0.050 

1  960 

3.510 

5.320 

8.510 

11 .710 

13.690 

STATIOI 

9 

x 

«  160.20 

lEDDY 

•  3 

Y- 

0.000 

0.490 

8.900 

14.660 

17.430 

19.360 

20.390 

20.828 

Z- 

0.025 

0.025 

1  .935 

3.515 

5.185 

7.555 

10.605 

13.690 

STATIOI 

10 

I 

•  178.00 

lEDDY 

«  3 

Y- 

0.000 

0.460 

12 .670 

16  510 

18.660 

19  730 

20.510 

20.899 

z- 

0.000 

0  000 

2.650 

4.160 

5.980 

7.690 

10.030 

13.690 
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SHIP  OFFSETS  (II  FT)  iT  71  STiTIOIS 
SCILIIG  FACTORS  :  Y  «>  1.000  2  «>  1.000 

TRIMHED  OFFSETS  FRQH  HIOSHIPS  DRAFT  AID  TRIH 


STATIOI 

11 

1 

•  195.80 

lEDDY  ■  3 

Y- 

0.000 

0.430 

12.670  16.390 

18.530 

19.540 

20.360 

20.815 

2- 

-0.025 

-0.025 

2.685  4.295 

6.135 

7.635 

10.125 

13.690 

STATIOI 

12 

I 

>  213.60 

lEDDY  ■  3 

Y- 

0.000 

0.430 

9.900  14.330 

17.400 

18.870 

20.080 

20.692 

2- 

-O.OSO 

-0.050 

2.090  3.600 

5.610 

7.230 

9.850 

13.690 

STATtOI 

13 

X 

•  231.40 

lEDDY  ■  3 

0.000 

0.430 

11.530  14.720 

16.690 

18.740 

19.680 

20.539 

2- 

-0.075 

-0.075 

2.915  4.265 

5.685 

7.825 

9.795 

13.690 

STATIOI 

14 

X 

•  249.20 

lEDDY  -  3 

Y- 

0.000 

0.400 

9.560  13.950 

17.480 

18.760 

19.680 

20.301 

2- 

0.120 

0.120 

2.720  4.640 

7.240 

8.950 

11 .040 

13.690 

STATIOI 

15 

X 

«  267.00 

lEDDY  •  2 

Y- 

0.000 

0.400 

8.360  14.560 

16.950 

18  510 

19 . 390 

19.741 

2- 

0.335 

0.835 

3.275  6.045 

7.865 

10.015 

12.285 

13.690 

STATIOI 

16 

X 

«  284.80 

lEDDY  «  2 

Y- 

0.000 

0.510 

8.410  11.210 

14.580 

16.450 

17.890 

18.697 

2- 

2.550 

2.550 

4.780  5.820 

7.510 

9.110 

11.200 

13.690 

STATIOI 

17 

X 

»  302.60 

lEDOY  ■  2 

Y- 

0.000 

0.460 

7.660  12.550 

14.650 

16.350 

17.020 

17.412 

2- 

5.235 

5.235 

6.885  8.375 

9.815 

11.385 

12.685 

13.690 

STATIOI 

18 

X 

•  320.40 

lEDDY  »  2 

Y- 

0.000 

0.380 

10.200  12.260 

13.550 

14.620 

15.380 

15.478 

2- 

8.260 

8.280 

9.860  10.530 

11.280 

12.220 

13.430 

13.690 

STATIOI 

19 

X 

»  338  20 

lEDDY  »  0 

Y- 

0.000 

6.560 

8.930  10.220 

11.330 

11.970 

12.610 

12.778 

2- 

10.435 

10.525 

10.665  11.145 

11 .795 

12  425 

13.365 

13.690 

STATIOI 

20 

X 

•  356.00 

lEDDY  •  0 

Y- 

0.030 

4.740 

6.250  7.230 

6.150 

8.800 

9  650 

9  794 

2- 

10.410 

10.560 

10.720  11.070 

11.600 

12.200 

13.380 

13.690 
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BILGE  lEEL  DATA  (II  FT) 


■0 

OF  BILGE  KEEL  PAIRS 

■  I 

BK 

STATIOI 

YBI 

ZBX 

BREADTH 

LEIGTH 

1 

10 

18.960 

7.010 

2.000 

12.100 

1 

11 

18.760 

6.100 

2.000 

17.800 

1 

12 

18.470 

6.000 

2.000 

17.800 

1 

13 

18.250 

6.380 

2.000 

17.800 

1 

14 

18.210 

7.130 

2.000 

21 .900 

*•••••••••*•««*•••  SIEG  DATA  (II  FT)  •••••*»•••••••»••*• 

STATIQI  «  17.87  BREADTH  -  3.900  LEIGTH  >60.000 


foil  data  (  FPS  UIITS  ) 

ID  OF  FII  PAIRS  «  0 

10  OF  SHAFT  BRACKET  PAIRS  ■  2 


FOIL 

STATIOI 

YF 

ZF 

SPAI 

RT  CHORD 

TP  CHORD 

CLA 

DI  AIGLE 

VEFF 

1 

>  BRACKET 

19.30 

10 . 200 

10.500 

7.400 

2.625 

2.625 

0.000 

79.000 

1.000 

2 

>  BRACKET 

19.30 

3.800 

9.700 

7.350 

2.625 

2.625 

0.000 

65.000 

1 .000 

3 

>  RUDDER 

19.16 

6.500 

5.850 

9.920 

7.500 

6.017 

0.000 

90.000 

1.000 

10  OF  RUDDERS  «  2  /  II  PROP  SLIPSTREAH  :  YES 


COITROLS 

ACCEL  GAII 

VEL  GAII 

GAII 

lATURAL  FREQ 

DAHPIIG  RATIO 

HIFREQ  CUTOFF 

LOFREQ  CUTOFF 

RUDDER  ROLL 

0.000 

0.000 

0.000 

0.000 

0.0000 

0.0000 

0.0000 

RUDDER  YAH 

0.000 

0.000 

0.000 

ROLL  RADIUS  OF  GYRATIOI  >  4.04  H 

RRG/BEAH  >  0.3169 

RRG/R  «  0.4793  WHERE  R«SqRT((B/2)«»2+T«»2) 
ROLL  lATURAL  FREQUEICY  «  0.560S 


59 


C.Q.  FftEQUEICY  RESPOISE  tl  REGULM  VtVES 
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ROLL  DIMPIIC  COEFFICIEnS  II  lEGtUR  HIVES 


DiltPIlG  COEFFICIEITS  101-DIItEISIQIU.IZED  BY  HUITIPLYIIG  BY  llIF/(2 . *044) 

(i.a.  roll  doiping  glron  aa  fraction  of  critical) 

ADDED  BOLL  IIEITII  lOI-DIKEISIOIALIZED  BY  ROLL  IIERTIA 

UF  ■  0.661  RAD/SEC  C44  ■  0.1833Et08  l-M  144  •  0.4778E+08  IG-II«»2 

ROUIG  IS  TIE  ROLL  AIGLE  II  DEGREES  AT  HIICI  TIE  ROLL  DAKPIIG  IAS  BEEI  COMPUTED 

SEA  DIRECTIOI  -  0.0  DEG 

U  •  18.0  UOTS  FROUDE  10  >  0.284 


SEA  STATE  >5,  SIG  HAVE  IT  «  3.249  M  ,  HAVE  PERIOD  •  10.000  SEC  ,  SPECTRUM  >  BRETSCHIEIDER 


HH 

HE 

ROUIG 

B44T0TiL 

BLIFT 

BHAVE 

BBI 

BEDDY 

BVBUU 

BVAPPEID 

A44T0TAL 

0.200 

0.162 

1.3 

0.137 

0.119 

•  0.000 

O.OIS 

0.000 

0.001 

0.000 

!  0.205 

0.300 

0.216 

1.3 

0.140 

0.119 

•  0.000 

0.020 

0.000 

0.001 

0.000 

!  0.206 

0.400 

0.249 

1.3 

0.141 

0.119 

!  0.000 

0.022 

0.000 

0.001 

0.000 

■  0.207 

0.600 

0.264 

1  .3 

0.141 

0.119 

!  0.000 

0.022 

0  000 

0.001 

0.000 

!  0.207 

0.600 

0.260 

1.3 

0.141 

0.119 

!  0.000 

0.022 

0.000 

0.001 

0.000 

!  0.207 

0.700 

0.237 

1.3 

0.140 

0.119 

•  0.000 

0.021 

0.000 

0.001 

0.000 

!  0.206 

0.800 

0.196 

1.3 

0.139 

0.119 

•  0.000 

0.019 

0.000 

0.001 

0.000 

!  0.206 

0.900 

0.135 

1.3 

0.136 

0.119 

•  0.000 

0.017 

0.000 

0.001 

0.000 

•  0.205 

1.000 

O.OSS 

1.3 

0.134 

0.119 

!  0.000 

0.015 

0.000 

0.001 

0.000 

!  0.204 

1.100 

0.043 

1 .3 

0.134 

0.119 

f  0.000 

0.015 

0.000 

0.002 

0.000 

•  0.204 

1.200 

0.161 

1.3 

0.137 

0.119 

!  0.000 

0.018 

0.000 

0.001 

0.000 

•  0.206 

1.300 

0.297 

1.3 

0.143 

0.119 

!  0.000 

0.023 

0.000 

0.001 

0.000 

!  0.208 

1.400 

0.452 

1.3 

0.148 

0.119 

•  0.001 

0.028 

0.000 

0.001 

0.000 

!  0.214 

1.600 

0.626 

1 .3 

0.156 

0.118 

!  0.006 

0.032 

0.000 

0.001 

0.000 

!  0.223 

1.600 

0.819 

1 .3 

0.174 

0.117 

!  0.020 

0.035 

0.000 

0.001 

0.000 

!  0.228 

1.700 

1 .031 

1.3 

0.204 

0.116 

•  0.047 

0.039 

0.000 

0.001 

0.001 

!  0.217 

1.800 

1.262 

1 .3 

0.232 

0.115 

•  0.072 

0.042 

0.000 

0.001 

0.001 

!  0.186 

1.900 

1.511 

1 .3 

0 . 2:  8 

0.114 

•  0.077 

0.046 

0.000 

0.001 

0.001 

!  0.164 

2.000 

1  .780 

1 .3 

0.230 

0.112 

•  0.066 

0.049 

0.000 

0.001 

0.002 

•  0.139 
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C  O.  FREQUaCT  lESPOISE  II  REGULil  HIVES 
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RuiX  DIHPIIC  COEFFICIEITS  II  lEGULAR  WIVES 

DIMPIIG  COEFFICIEITS  lOI-DDCElSIOIlLIZED  BY  NULTIPLYIIG  BY  RIF/(2.*C44) 

(i.«.  roll  daMpiog  gi**n  as  fraction  of  critical) 

ADDED  HOU.  IIERTIl  lOI-DIMElSIOIlLIZED  BY  ROLL  IIERTIl 

RIF  >  0.S61  RIO/SEC  C44  •  0.1832E-K>6  l-H  144  -  0.4TT6E*0e  IG-R**2 

ROLIIG  IS  THE  ROLL  ilGLE  II  DEGREES  IT  WHICH  THE  ROLL  DIMPIIG  HAS  BEEI  COMPUTED 

SEA  DIRECTIOI  •  30.0  DEG 

U  >  18.0  RIOTS  FROUDE  10  «  0.284 


SEA  STATE  -  5,  SIG  WAVE  HT  =  3.249  M  ,  WAVE  PERIOD  »  10.000  SEC  ,  SPECTRUM  *  BRETSCHIEIDER 
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C.G.  FXEQUEICY  HESPOISE  II  REGULM  HIVES 
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ROLL  DIMPIIG  COEFriCIEITS  II  REGULiR  HIVES 


DIHPIIG  COEFEICIEITS  lOI-DIKEISiafiLIZED  BY  HULTIPLYIIG  BY  UF/(2.*C44) 
(i.«.  roll  duipiiig  gi*«n  fraction  of  critical) 

ADDED  ROLL  IIERTIl  lOl-DIXEISIOBALIZED  BY  ROLL  IlERTIA 

RIF  «  0.S61  RID/SEC  C44  •  0.ia33E'»08  I-M  144  «  0.4778E+08 

ROLAIG  IS  THE  ROLL  AIGLE  II  DEGREES  AT  HBICE  THE  ROLL  DAHPIIG  BAS  BEEI  COKPUTED 

SEA  DIRECTIOI  -  60.0  DEG 

U  •  18.0  UQTS  FROUDE  10  >  0.284 


SEA  STATE  «  S,  SIG  HAVE  HT  ■  3.249  H  .  WAVE  PERIOD  «  10.000  SEC  ,  SPECTRUM  «  BRETSCHIEIDER 


HU 

HE 

ROLAIG 

B44T0TAL 

BLIFT 

BWAVE 

BBK 

BEDDY 

BVHULL 

BVAPPEID 

A44T0TAL 

0.200 

0.161 

!  10.2 

0.165 

0.119 
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0.001 
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0.001 
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C.G.  FREqUEICT  RESPOISE  II  REGULIR  URVES 
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ROLL  DIMPIIG  COEFFICIEITS  II  REGULAR  HAVES 


DAMPIIG  COEFFICIEITS  lOI-DIXEISIOIALIZED  BV  NULTIPLTIIC  BY  RIF/(2.*C44) 

(i.*.  roll  dtaping  glTon  as  fraction  of  critical) 

ADDED  ROLL  IIERTIA  lOI-DIXEISIOIALIZED  BY  ROLL  IIERTIA 

RIF  -  O.sei  RAD/SEC  C44  •  0.1832E«08  I-M  144  >  0.4778E'K>8  1G-M«*2 

ROUIG  IS  THE  ROLL  AIGLE  II  DEGREES  AT  HIICH  TIE  ROLL  DAHPIIG  HAS  BEEI  COHPUTED 

SEA  DIRECTIOI  >  90.0  DEG 

U  •  18.0  UOTS  FROUDE  10  >  0.284 


SEA  STATE  •  B,  SIG  HAVE  HT  •  3.249  H  .  HAVE  PERIOD  •  10.000  SEC  .  SPECTRUM  «  BRETSCHIEIDER 
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B44T0TAL 

BLIFT 
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BBX 

BEDDY 

BVHULL 

BVAPPEID 

A44T0TAL 

0.200 

0.200 

2.8 

0.147 

0.119 

0.000 

0.027 

0.000 

0.001 

0.000 

•  0.135 

0.300 

0.300 

2.8 

0.152 

0.119 

0.000 

0.032 

0.000 

0.001 

0.000 

!  0.135 
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2.8 
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0.037 
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!  0.135 
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0.054 
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!  0.135 
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0.056 
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!  0.135 
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0.243 
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0.000 
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0.002 

!  0.135 
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1.300 

2.8 

0.253 

0.115 

0.075 

0.060 

0.000 

0.001 

0.002 

!  0.135 
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1 .400 

2.8 
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0.001 
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0.003 

!  0.135 
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!  0.135 
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0.004 

!  0.135 
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2.8 
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0.001 

0.001 

0.004 

•  0.135 
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2.8 
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0.065 

0.068 

0.001 

0.001 
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!  0.135 
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0.001 

0.001 
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!  0.135 
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2.000 

2.8 

0.246 

0.111 
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0.001 

0.006 
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C.Q.  FKEqUEICY  RESPOISE  II  RECULM  WIVES 
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4t  J«H(^$«>4C<4C0iA<rc«T4OOOOO5oOO 

a,  •^^ddv^dddddddddddddd 
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(A«^««^M«PP«9c<MnkApf4nNPP 

wcdddddddddddNsd^eoNd 

IdCL  lllill*«kAKP  nvCiVCMpCN 
>  I  I  I  I  ^  ^  ^  I  ^ 


r)r««rtA«Kr<»oc«9i«*kAppr)MVKfO 

8PPPPPC^P«-«v^PKkA^v^v>^v«PO 

ppff>ppp»^eo«ppoppppoo 


»<*-i-OOP^*-PPPPPOPPPPP 


ba9>p«opnpp^pnr«^0>p^pr(c«p 

vj^o^e*®pp*<KOK®^OP^nkAO 

^ . 

CB®®^NPOPvr«P®N«kA®^PPN 
0>®®0>®®®®POC4iA<^®r>  I  00®h<. 


NPP®lA^r)C<^OPOOPOOPOO 

ddddddddddddddddddo 


u^v^K®c4p«^piAv4«pr««4eDno^h> 

MOON^tAPKVO^nnOVNOPDh- 

. . 

B®<r®iA<*)(^P9>®®r«®N®PS>9>^P 

wa.«^OP®o>pp®cD®®®o®®oiAvp 

®  ^  ^  I  I  I  ^  ^  ^ 

at  I 

p 

V}  kA®^r<o>f*)®^r)iAr«ppo®r}c<c<r4 
^^iAr«®tAr«®^0®r)pptMppoop 
E^r)r>r*e*r«»-»^*i4ppppppp555 

ddddddddddddodddddd 


pno®oc«r«f*)r«e«p®p(OOtAv^pp 

v-'^r«.^Kno®KK»pc<o^®r)OP 
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ROLL  DIKPIIG  COEFFICIEITS  II  RECULIR  VIVES 

DIKPIIG  COEFFICIEITS  IOI-DIIIEISIQIU.IZED  BY  HULTIPLYIIC  BY  UF/(2.*C44) 

(i.«.  roll  dOBpiag  giros  ••  fraction  of  critical) 

ADDED  ROLL  IIERTIA  IQI-DDIEISIOIALIZED  BY  ROU.  IIERTll 

UF  •  0.561  UD/SEC  C44  -  0.1833E«08  B~H  144  •  0.4778E'K>e  IG-H**2 

ROLIIG  IS  TIE  ROU.  AIGLE  II  DEGREES  AT  VBICB  TIE  ROU.  DAIIPIIG  BAS  BEEI  COMPUTED 

SEA  DIRECTIOI  -  120.0  DEC 

U  •  18.0  UOTS  FROUDE  10  «  0.284 


SEA  STATE  •  S.  SIG  HAVE  BT  >  3.249  H  ,  HAVE  PERIOD  •  10.000  SEC  ,  SPECTRUM  >  BRETSCBIEIDER 


w 

HE 

ROUIG 

B44T0TAL 

BLIFT 

BVAVE 

BBK 

BEDDY 

BVBUU 

BVAPPEID 

A44T0TAL 

0.200 

0.219 

2.3 

0.146 

0.119 

0.000 

0.036 

0.000 

0.001 

0.000 

!  0.143 

0.300 

0.343 

2.3 

0.151 

0.119 

0.000 

0.032 

0.000 

0.001 

0.000 

!  0.143 

0.400 

0.476 

3.3 

0.157 

0.118 

0.001 

0.037 

0.000 

0.001 

0.000 

!  0.143 

0.500 

0.618 

2.3 

0.165 

0.118 

0.005 

0.041 

0.000 

0.001 

0.000 

•  0.143 

0.600 

0.770 

2.3 

0.179 

0.118 

0.016 

0.045 

0.000 

0.001 

0.000 

!  0.143 

0.700 

0.932 

2.3 

0.200 

0.117 

0.033 

0.048 

0.000 

0.001 

0.001 

•  0.143 

0.800 

1.102 

2.3 

0.226 

0.116 

0.056 

0.053 

0.000 

0.001 

0.001 

!  0.143 

0.900 

1.283 

2.3 

0.247 

0.115 

0.074 

0.055 

0.000 

0.001 

0.002 

•  0.143 

1.000 

1 .472 

2.3 

0.254 

0.114 

0.077 

0.059 

0.000 

0.001 

0.003 

!  0.143 

1.100 

1.672 

2.3 

0.251 

0.113 

0.071 

0.062 

0.001 

0.001 

0.003 

!  0.143 

1.200 

1.880 

2.3 

0.244 

0.112 

0.062 

0.065 

0.001 

0.001 

0.004 

!  0.143 

1.300 

2.098 

2.3 

0.238 

0.110 

0.052 

0.068 

0.001 

0.001 

0.006 

;  0.143 

1.400 

2.336 

2.3 

0.232 

0.109 

0.044 

0.071 

0.001 

0.001 

0.007 

!  0.143 

1.500 

2.563 

2.3 

0.225 

0.107 

0.034 

0.074 

0.002 

0.001 

0.008 

!  0.143 

1.600 

2.810 

2  3 

0.334 

0.106 

0.029 

0.077 

0.002 

0.001 

0.010 

!  0.143 

1.700 

3.065 

2.3 

0.320 

0.104 

0.032 

0.079 

0.002 

0.001 

0.011 

!  0.143 

1.800 

3.331 

2  3 

0.222 

0.103 

0.021 

0.082 

0.002 

0.001 

0.013 

!  0.143 

1.900 

3.606 

2.3 

0.225 

0.101 

0.021 

0.085 

0.003 

0.001 

0.014 

•  0.143 

2.000 

3.890 

2.3 

0.318 

0.100 

0.011 

0.088 

0.003 

0.001 

0.016 

<  0.143 
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C.6.  FREqUEICT  KESPOISE  II  REGULIR  HiVES 
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EOLL  DiHPIlG  COEFFICIEITS  II  lEGUUR  HAVES 


DAHPIIG  COEFFICIEITS  lOI-OIXEISIOIALIZED  BY  HULTIPLYIIC  BY  MF/(2.*C44) 

(i.«.  roll  (taping  given  u  fraction  of  critical) 

ADDED  lOU.  IIEBTIA  lOI-DDIEISIOIAlIZED  BY  lOU.  IIERTIA 

UF  -  0.561  UD/SEC  C44  ■  0.1832E*08  l-R  144  >  0.477BK*0B  1C-H**3 

lOUIG  IS  TIE  BOU  AIGLE  II  DEGREES  AT  VHICI  TIE  RQU.  DAHPIIG  IAS  BEEI  COHPIH'EO 

SEA  DIRECTIOI  •  ISO.O  DEG 

U  ■  18.0  UOTS  FROUDE  10  •  0.384 


SEA  STATE  >  5,  SIG  HAVE  IT  •  3.249  M  ,  HAVE  PERIOD  ■  10.000  SEC  ,  SPECTRUM  •  BRETSCHIEIDER 


HH 

HE 

ROLAIG 

B44T0TAL 

BLIFT 

BHAVE 

BBX 

BEDDY 

BVEUIX 

BVAPPEID 

A44T0TAI. 

0.200 

0.233 

1 .3 

0.141 

0.119 

0.000 

0.021 

0.000 

0.001 

0.000 

•  0.150 

0.300 

0.374 

1 .3 

0.145 

0.119 

0.000 

0.026 

0.000 

0.001 

0.000 

•  0.150 

0.400 

0.531 

1  .3 

0.151 

0.118 

0.003 

0.030 

0.000 

0.001 

0.000 

!  0.150 

0.500 

0.705 

1.3 

0.163 

0.118 

0.010 

0.034 

0.000 

0.001 

0.000 

!  0.150 

0.600 

0.695 

1  .3 

0.183 

0.117 

0.028 

0.037 

0.000 

0.001 

0.000 

!  0.150 

0.700 

1.101 

1 .3 

0.214 

0.116 

0.055 

0.040 

0.000 

0.001 

0.001 

•  0.150 

0.800 

1  .324 

1.3 

0.236 

0.115 

0.075 

0.044 

0.000 

0.001 

0.001 

'  0.150 

0.900 

1.563 

1.3 

0.238 

0.114 

0.075 

0.047 

0.000 

0.001 

0.002 

!  O.ISO 

1.000 

1.818 

1.3 

0.230 

0.112 

0.064 

0.050 

0.000 

0.001 

0.002 

!  0.150 

1.100 

2.090 

1.3 

0.220 

0.110 

0.052 

0.053 

0.001 

0.001 

0.003 

!  O.ISO 

1.200 

2.378 

1.3 

0.211 

0.108 

0.042 

0.056 

0.001 

0.001 

0.004 

!  0.150 

1.300 

3.683 

1 .3 

0.202 

0.106 

0.030 

0.059 

0.001 

0.001 

0.005 

•  0.150 

1.400 

3  004 

1 .3 

0.198 

0.104 

0.024 

0.062 

0.001 

0.001 

0.006 

!  0.150 

1.500 

3.341 

1 .3 

0.197 

0.103 

0.021 

0.064 

0.001 

0.001 

0.007 

!  0.150 

1.600 

3.695 

1.3 

0.199 

0.101 

0.020 

0.067 

0.002 

0.001 

0.008 

!  0.150 

1.700 

4.065 

1 .3 

0.193 

0.099 

0.012 

0.070 

0.002 

0.001 

0.009 

•  0.150 

1.800 

4.451 

1  .3 

0.190 

0.097 

0.006 

0.072 

0.002 

0.001 

0.011 

•  0.150 

1.900 

4  854 

1  .3 

0.193 

0.096 

0.006 

0.075 

0.002 

0.001 

0.012 

!  0.150 

2.000 

5.273 

1  .3 

0.197 

0.094 

0.008 

0.078 

0.003 

0.001 

0.013 

•  0.150 
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ROLL  DiHPIlG  COEFFICIEITS  II  REGULAR  HAVES 


DAHPIIG  COEFFICIEITS  lOI'OIXEISIOIALIZEO  BY  HULTIPLYIIC  BY  RIF/(3.*C44) 

(i.«.  roll  daaping  glTOD  u  fraction  of  critical) 

ADDED  ROU  IIERTIA  lOI-DIHEISIOIALlZED  BY  ROLL  IIERTIA 

UF  >  O.S61  RAD/SEC  C44  ■  O-lBSlE-KiB  l-H  144  -  0.4778E+08  IC-H*«2 

ROUIQ  IS  TIE  ROLL  AIGLE  II  DEGREES  AT  VBICB  TIE  ROLL  DAHPIIG  IAS  BEEI  COHPUTED 

SEA  DIRECTIOI  >  180.0  DEG 

U  ■  18.0  UOTS  FROUDE  10  -  0.384 


SEA  STATE  >  5,  SIG  HAVE  HT  <•  3.249  H  .  HAVE  PERIOD  >  10.000  SEC  ,  SPECTRUM  «  BRETSCBIEIDER 


HH 

HE 

ROUIG 

B44T0TAL 

BLIFT 

BHAVE 

BBK 

BEDDY 

BVaULL 

BVAPPEID 

A44T0TAL 

0.200 

0.338 

1.3 

0.140 

0.119 

0.000 

0.021 

0.000 

0.001 

0.000 

!  0.206 

0.300 

0.385 

1 .3 

0.145 

0.119 

0.000 

0.026 

0.000 

0.001 

0.000 

•  0.212 

0.400 

0.551 

1  .3 

0.152 

0.118 

0.003 

0.030 

0.000 

0.001 

0.000 

!  0.220 

0.500 

0.736 

1.3 

0.165 

0.118 

0.013 

0.034 

0.000 

0.001 

0.000 

.'  0.227 

0.600 

0.940 

1 .3 

0.190 

0.117 

0.035 

0.037 

0.000 

0.001 

0.000 

!  0.224 

0.700 

1.163 

1 .3 

0.221 

0.116 

0.063 

0.041 

0.000 

0.001 

0.001 

!  0.200 

0.800 

1.405 

1  .3 

0.238 

0.114 

0.077 

0.044 

0.000 

0.001 

0.001 

!  0.166 

0.900 

1.665 

1 .3 

0.235 

0.113 

0.071 

0.048 

0.000 

0.001 

0.002 

!  0.144 

1.000 

1 .945 

1 .3 

0.224 

0.111 

0.059 

0.051 

0.000 

0.001 

0.003 

•  0.136 

1.100 

3.343 

1.3 

0.215 

0.109 

0.047 

0.054 

0.001 

0.001 

0.003 

!  0.134 

1.200 

3.561 

1.3 

0.204 

0.107 

0.034 

0.057 

0.001 

0.001 

0.004 

!  0.136 

1.300 

3.897 

1.3 

0.200 

0.105 

0.028 

0.060 

0.001 

0.001 

0.005 

•  0.139 

1.400 

3.353 

1.3 

0,197 

0.103 

0.022 

0.063 

0.001 

0.001 

0,007 

!  0.142 

1.500 

3.636 

1.3 

0.198 

0.101 

0.021 

0.066 

0.001 

0.001 

0.008 

•  0.146 

1.600 

4.019 

1 .3 

0.192 

0.099 

0.013 

0.069 

0.002 

0.001 

0.009 

•  0.145 

1.700 

4.431 

1 .3 

0.188 

0,097 

0.005 

0.071 

0.002 

0.001 

0.010 

!  0.145 

1.800 

4.862 

1 .3 

0.191 

0,096 

0.006 

0.074 

0.002 

0.001 

0.012 

!  0.148 

1 .900 

5.311 

1.3 

0.195 

0,094 

0.008 

0.077 

0.002 

0.001 

0.013 

!  0.150 

2.000 

5,780 

1  .3 

0.194 

0.092 

0.0C3 

0.080 

0.003 

0.001 

0.015 

!  0.151 
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FORCES  RELiTIVE  TO  LOCIL  IXES  110  NOTIOI-IIDUCED  IRTERRUPnOIS 

LIT  TIP  COEFF  *0.250  ,  LOIG  TIP  COEFF  -  0.170  ,  TIKE  FOR  OPERITIDI  •  60.0  s 


BEiDIBG 

••*LiTERlL*** 

**LQIGITUDIIAL** 

TOTAL 

LFE 

Nil 

LFE 

Mil 

Nil 

DEG 

G 

G 

0.0 

0.000 

0.000 

0.001 

0.000 

0.000 

30.0 

0.061 

0.001 

0.002 

0.000 

0.001 

60.0 

0.139 

1.938 

0.004 

0.000 

1.938 

90.0 

0.077 

0.150 

0.001 

0.000 

0.150 

130.0 

0.063 

0.138 

0.025 

0.130 

0.267 

ISO  0 

0.028 

0.003 

0.028 

0.214 

0.217 

180.0 

0.000 

0.000 

0.028 

0.189 

0.189 

sumiiG 

UD  DECK  HETIESS  CILCUUTIOIS  - 

'  FREEBOARD 

AT  STATIOI 

FREEBOiRO 

AFTER  HAVE 

PROFILE 

CORRECTIOl 

IS  S.7S 

H 

lEADIlG 

DEG 

••••VERTICAL****** 

REL  DISP  REL  VEL 

H  H  /SEC 

VRgi 

PRQB 

(DV) 

DU  PER 

BOOR 

PROB 

<KE) 

IE  PER 

BOUR 

SLAN  PRESSURE 
HOSTPROB  EXTREME 

IPA  IPA 

SLAN  FORCE 
HOSTPROB  EXTREME 

I  /  M  •  /  M 

0.0 

0.748 

0.215 

0.003 

0.0000 

0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

30.0 

0.712 

0.203 

0.004 

0.0000 

0.0 

0,0000 

0.0 

0.00 

0.00 

0.0 

0.0 

60.0 

0.559 

0.281 

0.016 

0.0000 

0.0 

0,0000 

0.0 

0.00 

O.oO 

0.0 

0.0 

90.0 

0.154 

0.245 

0.076 

0.0000 

0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

120.0 

1  .129 

1.587 

0.211 

0.0000 

0.0 

0,0013 

1.0 

37.41 

94.25 

33628 . 9 

85221.0 

150.0 

1.427 

1 .984 

0.217 

0.0003 

0.2 

0.0155 

12.3 

106.41 

195.29 

96217.4 

176578.5 

180.0 

1.475 

2.C78 

0.212 

0.0005 

0.4 

0.0203 

16.3 

122.68 

220.14 

110925 .0 

199044.9 
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MS  ROTIOIS  II  SROKT-CKESTED  SEAS  -  SPRElDIia  IIGLE  *  90.00  DEQ 
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FORCES  RaiTIVE  TO  LOCAL  IIES  ilD  BOTIOI-IIDUCFO  IITE8RUPTIDIS 

LiT  TIP  COEFF  -  0.250  .  LOIG  TIP  COEFF  -0.170  ,  TME  FOR  OPERITIOI  -  60.0  » 


BEADIIG 

••-LATERAL*** 

••LQIGITUDIIAL** 

TOTAL 

LFE 

mi 

LFE 

HII 

mi 

DEG 

G 

G 

0.0 

0.071 

0.017 

0.003 

0.000 

0.017 

30.0 

0.082 

0.100 

0.003 

0.000 

0.100 

60.0 

0.095 

0.411 

0.007 

0.000 

0.411 

90.0 

0.090 

0.427 

0.015 

0.000 

0.427 

120.0 

0.068 

0.104 

0.022 

0.020 

0.124 

150.0 

0.043 

0.008 

0.026 

0.123 

0.131 

180.0 

0.033 

0.003 

0.028 

0.191 

0.193 

sumiiG 

AID  DECK  VETIESS  CALCimATIOIS  - 

FREEBOARD 

AT  STATIOI  IS 

6.03  H 

FREEBOIRD  AFTER  VIVE 

;  PROFILE 

CORRECTIOI 

IS  5.75 

H 

BEADIIG 

••••VERTICAL****** 

VRQI  PRQB  DV  PER 

PROB 

IE  PER 

SLAB  PRESSURE 

SLAM  FORCE 

REL  DISP  REL  VEL 

(OV)  HOUR 

(IE) 

HOUR 

MOSTPROB 

EXTREME 

MOSTPROB 

EXTREME 

DEG 

H  R 

[  /SEC 

EPA 

IPA 

I  /  H 

I  /  M 

0  0 

0.703 

0.222 

0.008  0.0000  0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

30.0 

0.658 

0.242 

0.025  0.0000  0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

60.0 

0.616 

0.506 

0.072  0.0000  0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

90.0 

0.773 

0.990 

0.130  0.0000  0.0 

0.0000 

0.0 

0.00 

0.00 

0.0 

0.0 

120  0 

1  .064 

1.478 

0.178  0.0000  0  0 

0.0006 

0.4 

23.33 

72.65 

21094.3 

65684 . 1 

150  0 

1.307 

1.828 

0.206  0.0001  0.1 

0.006S 

5.6 

77.32 

152.79 

69910.1 

138152.3 

180.0 

1  .400 

1 .958 

0  214  0.0002  0.2 

0.0132 

10.6 

100.74 

187.29 

91084.1 

169342.0 
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